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ABSTRACT 


\, 

^Accurate  ocean  data  is  essential  for  successful  fleet  operation.  The  N-ROSS 
Satellite,  which  is  being  developed  for  this  mission,  will  carry  a  Low  Frequency 
Microwave  Radiometer  (  LFMR  ).  The  LFMR  consists  of  large  flexible  reflector  and 
boom  and  spins  at  15  r.p.m.  The  effects  of  the  flexibility  of  the  boom,  the  spin-up 
procedure  and  the  structural  damping  on  the  pointing  error  of  the  LFMR  are 
investigated  by  performing  the  dynamic  simulation  using  the  Dynamic  Simulation 
Language.  Two  cases  of  boom  material,  Aluminum  Alloy  and  the  Graphite/epoxy 
composite  material,  are  analyzed  and  the  results  are  compared.  The  simulation  and 
analysis  results  are  presented  in  graphical  forms. 
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I.  INTRODUCTION 


A.  BACKGROUND 

Accurate  ocean  weather  prediction  is  essential  for  successful  fleet  operations.  The 
NANY  needs  superior  data  collection  capability  to  obtain  the  data  density  and 
reliability  necessary  to  produce  consistently  accurate  forecasts  and  oceanographic  data. 
However  present  prediction  models  are  limited  by  the  quality  and  quantity  of  input 
data  which  come  mainly  from  ships.  The  obvious  approach  to  satisfy  these  purpose 
can  be  derived  from  satellite  observations.  Therefore,  the  NAVY  planned  the 
construction  of  the  Navy  Remote  Ocean  Sensing  System  (  N  -  ROSS  ).  [Ref.  I]  This 
system  consists  of  satellites  which  scan  the  earth  surface  and  provide  the  fleet  with 
timely  worldwide  knowledge  of  ocean  data  such  as  seasurface  wind  speed,  wind 
direction,  seasurface  temperature,  ice  edge  detection,  ocean  wave  height  and  ocean 
photograpy.  [Ref.  2[  To  satisfy  the  mission  requirements  the  N  -  ROSS  Satellite  will 
carry  several  sensors.  Among  these  sensors  the  Low  Frequency  Microwave  Radiometer 
(  LFMR  )  is  the  most  important  and  the  most  interesting  from  the  dynamics  of  the 
spacecraft  point  of  view.  The  function  of  the  LFMR  is  scanning  the  earth  surface  and 
measure  the  seasurface  temperature.  To  increase  the  scanning  area  the  deployable 
reflector  spins  at  15  r.p.m.  The  sizes  of  this  LFMR  reflector  and  boom  are  relatively 
large  compared  to  the  N  -  ROSS  Satellite  itself.  So  the  weight  of  this  boom  should  be 
light,  which  makes  the  boom  flexible.  By  this  reason,  there  exist  certain  extent  of 
deflection  at  the  tip  of  the  LFMR  boom  and  this  boom  vibrates  when  this  boom  is 
spinning.  Deflection  and  vibration  due  to  this  elastic  deformation  induce  pointing 
error  of  the  reflector  in  elevation  and  azimuth  angle.  However  there  is  strict  pointing 
error  requirement  of  the  LFMR.  [Ref.  2]  Therefore  analysis  of  the  flexible  LFMR 
boom  which  supports  the  sensor  payload  is  imperative  for  this  research. 

B.  STATEMENT  OF  PROBLEM 

The  traditional  approach  to  dynamics  of  a  boom  system  is  based  on  the 
assumption  that  the  systems  are  composed  of  rigid  bodies.  Until  recently,  only  a  rigid 
body  motion  was  assumed  for  the  analysis.  However,  the  flexible  system  includes  a 
small  elastic  deformation  as  well  as  a  large  motion.  These  small  elastic  deformations 
include  bending,  twisting  inu  axial  extension.  Development  of  a  dynamic  modei 


including  flexibility  demands  more  accuracy  for  the  system  responses.  Without 
considering  these  small  motion  of  a  boom,  we  cannot  expect  a  certain  accuracy  to 
maintain  a  spacecraft  attitude  and  pointing  control.  [Ref.  3]  Recently,  efforts  have 
been  made  to  control  maneuvers  of  mechanical  systems  which  can  not  be  adequately 
modeled  using  a  rigid  body  assumption  for  all  or  some  of  the  system  components, 
especially  in  the  fields  of  satellites,  [Ref.  4:  pp.  257-264) 

Therefore,  the  development  of  a  good  dynamic  model  of  a  flexible  system,  an 
efficient  dynamic  equations  formulation  method  and  a  good  dynamic  simulation  schem 
are  essential  for  the  analysis  of  and  identification  of  potential  problems  in  the  flexible 
LFMR  system. 

C.  THESIS  OUTLINE 

In  Chapter  II,  the  development  of  an  analytic  model  for  the  Lower  Frequency 
Microwave  Radiometer  (  LFMR  )  reflector  boom  in  3-dimensional  motion  is  described. 
The  large  motion  due  to  rotation  is  described  by  an  equivalent  rigid  boom  motion  and 
elastic  deflection  of  a  flexible  boom  relative  to  the  equivalent  rigid  boom  motion  is 
expressed  using  the  mode  superposition  technique.  The  dynamical  equations  for  this 
model  are  formulated  using  the  Lagrange's  method. 

In  Chapter  III,  The  computer  implementations  for  the  solution  of  the  obtained 
equations  are  explained.  For  the  modal  analysis  of  the  system.  NAsa  STRuctural 
ANalysis  (  NASTRAN  )  computer  program  was  used  and  Dynamic  Simulation 
Language  (  DSL  )  was  applied  to  solve  the  simultaneous,  nonlinear,  ordinary 
differential  equations.  The  LINPACK  subroutines  DGEFA  and  DGESL  are  also  used 
in  the  dynamic  simulation. 

In  Chapter  IV,  simulation  results  are  presented  to  investigate  the  deflection  and 
pointing  error  of  the  LFMR.  Compansons  are  made  by  changing  the  torque  input 
condition.  The  problems  considered  are  1)  the  effects  of  spin-up  procedure  on  the 
pointing  error  of  the  LFMR  reflector;  2)  the  effect  of  damping  on  the  settling  time  of 
pointing  error;  3)  the  equilibrium  configuration  of  LFMR  booms  due  to  constant 
rotating  speeds  For  the  comparison  purpose,  two  kinds  of  material,  aluminum  alloy 
and  composite  material,  were  assumed  as  the  LFMR  boom  materials. 

In  chapter  V,  Conclusions  are  made  from  the  research  and  some 
recommendations  for  future  work  in  the  area  of  the  dynamic  analysis  of  the  flexible 
LFMR  reflector  boom  svstem  are  iziven. 
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II.  FORMULATIONS  OF  DYNAMIC  EQUATIONS 
A  INTRODUCTION 

In  this  Chapter,  the  dynamical  equations  of  motion  for  the  LFMR  system  which 
rotates  in  three  dimensional  space  is  developed.  The  2-dimensional  planar  motion  of 
the  same  boom  was  also  studies  and  are  presented  in  Appendix  A.  A  simple  dynamic 
model  of  the  LFMR  system  is  developed  for  this  analysis.  The  Lagrangian  approach 
and  the  mode  superposition  technique  are  used  for  the  formulation  of  the  dynamic 
equations  of  the  flexible  LFMR  system. 

B.  DESCRIPTIONS  OF  THE  MODEL 

The  LFMR  shown  in  Fig.  2.1  consists  of  four  structures:  a  reflector,  an  upper 
reflector  boom,  a  lower  reflector  boom  and  an  electronics  box.  The  reflector  is  attached 
to  the  top  of  the  upper  reflector  boom.  The  upper  reflector  boom  and  the  lower 
reflector  boom  is  connected  by  a  boom  hinge.  The  electronics  box  is  attached  to  the 
bottom  end  of  the  lower  boom. 

For  our  analysis,  the  deployable  reflector  is  modeled  as  a  concentrated  mass  at 
the  tip  of  the  upper  reflector.  We  assume  the  boom  hinge  which  connects  the  two 
booms  is  stiff  and  firm  and  there  is  no  relative  motion  between  the  booms  after  the 
deployment  of  the  LFMR  boom.  Therefore  we  consider  the  whole  system  (  reflector, 
booms,  boom  hinge  )  as  a  one  body  system. 

The  LFMR  system  is  connected  to  the  Main  Bus  of  the  N  -  ROSS  Satellite  by  a 
Spacecraft  Boom.  The  attitude  of  the  N  -  ROSS  Satellite  is  controilec  by  a  Attitude 
Determination  And  Control  System  (  ADACS  )  very  accurately.  [Ref.  5]  Therefore,  it 
is  assumed  that  the  spin  axis  and  the  base  is  remain  fixed  in  tiie  reference  frame  fixed 
to  the  N  -  ROSS  Satellite.  The  N  -  ROSS  Spacecraft  moves  on  a  circular  orbit  with 
the  spin  axis  always  pointing  the  earth  center.  Hence,  the  gravitational  force  is  in 
equilibrium  with  the  centrifugal  force  in  the  orbit  plane:  the  LFMR  system  is  in  zero-g 
environment.  Therefore,  in  the  dynamic  model  of  the  present  studies,  the  reference 
frame  fixed  to  the  N  -  ROSS  Spacecraft  is  assumed  the  Newtonian  (  inertial  )  reference 
frame  and  the  LFMR  system  is  in  zero-g  environment. 

From  the  above  assumptions  the  dynamic  model  of  the  LFMR  system  is  defined 
as  shown  in  Figure  2.2.  The  giobai  coordinates  X.  V.  Z  is  fixed  in  the  inertial  reference 
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frame  and  a  moving  coordinate  x,  y,  z  (  local  coordinates  ).  which  is  a  body  Fixed 
coordinate  system,  is  defined  as  shown  in  Figure  2.2.  The  body  fixed  coordinate 
system  (  local  coordinates  )  is  attached  to  the  base  O.  The  local  z-axes  and  the  global 
Z-axes  are  the  common  axis  of  rotation.  The  local  x,  y-axes  and  the  global  X,  Y-axes 
are  in  the  same  plane  with  angle  difference  0. 

C.  LAGRANGE'S  EQUATION 

For  any  system  there  must  be  same  numbers  of  independent  coordinates  as  the  degrees 
of  the  freedom  of  the  system  to  completely  describe  the  motion  of  the  system.  The 
choice  of  coordinates  is  important  in  dynamic  analysis.  Such  independent  coordinates 
are  called  generalized  coordinates  and  are  denoted  by  the  letter  qf.  For  a  system  with  a 

set  of  n  independent  generalized  coordinates  qr  (  r  =  l.  2,  3, . n  ),  Lagrange's 

equations  are  expressed  as  [Ref.  &7] 


£_[iL  dT  +  d\J 

dt  dqk  <3qk  <3qk 


(eqn  2.1) 


(  k  =  1,  2.  3 . n  ) 

where  T  is  the  kinetic  energy,  U  is  the  potential  energy  and  is  the  generalized  forces 
which  is  defined  as  follows 


(eqn  2.2) 


The  dot  over  a  variable  means  derivative  with  respect  to  time.  F.  is  the  force  acting  on 
particle  j  and  Rj  is  the  instantaneous  position  of  particle  j  and  may  be  expressed  in 
terms  of  generalized  coordinates 


Rj  =  Rj  (  ^1’  ^2’  ^3* . %  }  (eqn 

and  6Rj  is  the  virtual  displacement  of  the  particle. 

To  describe  the  motion  of  the  LFMR  boom  system  which  composed  of  a  large 
slow  motion  due  to  rotation  and  a  small  fast  motion  due  to  elastic  vibration,  two  kinds 
of  generalized  coordinates  are  defined.  One  :s  9  for  rotation  of  boom  and  others  arc  qh 
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Figure  2.1  N-ROSS  Baseline  Configuration. 


Figure  2.2  Modeling  of  LFMR  boom  system. 

(  h  =  1,  2,  3 . .  n  )  for  h-th  mode  generalized  displacement,  where  n  is  number  of 

modes. 

Now  Lagrange's  equation  2.1  is  rewritten  as 


d  <JT  <3T 

"  io 


(eqn  2.4) 


and 


d  ,  dl  31  d\j  ^ 

——  l  -r r*  1  ”  -7—  +  -r —  "  Qh 
dt  d%  d%  d% 


(eqn  2.5) 


(  h  -  1,  2,  3, 


.  n ) 


The  external  force  acting  on  the  LFMR  system  is  assumed  the  torque  t  by  a 
torque  motor  at  the  root  of  the  boom.  The  contribution  of  this  torque  to  the 
generalized  forces  is  Qr  *  r  and  this  torque  does  not  contribute  to  Qjr  Damping 
forces  are  assumed  equal  to  the  modal  damping  value  in  this  analysis.  Since  the  modal 
damping  values  can  be  measured  easily.  Therefore  the  contribution  of  damping  forces 
to  the  generalized  forces  can  be  obtained  using  a  dissipation  function. 


D  *  _l_V2(;iaiiMi  q^t) 


(eqn  2.6) 


with 


Qh“  - 


3D 


*h 

-  2  5h  <nh  Mh  qh 


where 

^h: 

°v 


M 


h- 


modal  damping  ratio  of  h-rh  mode 
the  natural  frequency  of  h-th  mode 
the  modal  mass  of  h-th  mode 


Then  equation  2.4  and  2.5  finally  written  as 


d 

~d7 


31 

{w 


1  " 


dT 


<39 


<3U 

<30 


*  r 


(eqn  2.7) 


and 


d 

"dt" 


3  1 


<3T 

d% 


3\: 

d% 


XT  1  “XT  +  X-  ’  2  ;h  c,JhVIh  Mh 


iecin  2. Si 
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(  h  -  1,  2.  3 . n  ) 

D.  POSITION  AND  VELOCITY 

During  the  rotational  motion  of  the  LFMR  system,  the  boom  deforms.  The 
deformed  positions  of  a  generic  point  in  the  system  can  be  expressed  the  vector  sum  of 
a  vector  Rq(x)  from  the  origin  to  the  undeformed  position  of  the  point  and  a  vector 
W(x,t)  which  is  deflection,  as  shown  in  Figure  2.3.  The  notations  show  in  the  Figure 
2.3  represent  the  following  parameters: 
tip  mass 

mass  of  electronics  box 
length  of  lower  boom 
length  of  upper  boom 
angle  between  two  links 
applied  torque 
angular  displacement 
angular  velocity 

position  vector  of  the  point  on  the  boom  in  the  local  x-direction 
position  vector  of  the  point  on  the  boom  after  deformation 
deformation  vector  of  boom 
unit  vector  of  local  x-direction 
unit  vector  of  local  v-direction 
unit  vector  of  local  z-direction 
unit  vector  of  global  X-direction 
unit  vector  of  global  Y-direction 
unit  vector  of  global  Z-direction 

Then  the  position  vector  R(x.t)  is  expressed  as 


M  : 

mr; 

<1  = 

V  . 
P 

T  ; 

0(t>  : 
0(t)  : 

Rq(x>  : 
R(x.t)  : 
W(x.t)  : 
i  : 

j  : 
k  : 

‘o  : 

Jo ; 

N) : 


Rfx.t)  =  Rq(x)  -  W(x.t)  (eqn  2.9) 

The  undeformed  position  Rq(x )is  represented  by  its  components. 

R,,(x»  =  Rx(x)  +  R,(x)  ^eqn  2.10i 

=  Rx(.x)  i  +  Rz(x>  k 
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Figure  2.3  Parameters  of  the  boom  system. 

The  elastic  deflection  W(x,t)  is  expressed  as  the  modal  sum  as  follows: 

W(x.t)  -  V  [  <p(x(x)  i  +  <pi>'(x)  j  +  <p/(x)  k  ]  q.(t) 
i 


(eqn  2.1 1) 


where 


q>.*lx)  i  is  i-tk  mode  shape  function  in  extension 
9jy(x)  j  is  i-th  mode  shape  function  in  translation 

From  the  equation  2.10  and  equation  2.11,  equation  2.9  can  be  rewritten  in  the  form  of 


R(x,t)  *  Rx(x)  i  +  Rz(x)  k  (eqn  2.12) 

+  y  [  <P  X(x)  i  +  9f(x)  j  +  <p;z(x)  k  ]  q;(t) 
i 

=  [  Rx(x)  +  <px(x)  q/t)  ]  i  +  [  V  <py(x)  q;(t)  ]  j 

i 

+  [  Rz(x)  +  X  <PiZ(x)  q.(t)  ]  k 


The  velocity  of  the  point  in  the  Newtonian  reference  frame  is  obtained  by  taking  the 
time  derivative  of  equation  2.12  and  applying  the  relation  between  the  time  derivative 
of  unit  vector  i  and  j, 


i  -  e  k  x  i  =  e  j 

j  =  0  k  *  j  =  -  0  i 

it  a  0  k  x  k  ■  0 


the  velocity  is 

R(x.t)  =  [  -  B  y  <P;y(x)  q;( t )  +  y  <p.x(x)  q;(t)  ]  i  (eqn  2.13) 

l  l 

+  [  0  Rx(x)+  0  y  (PiX(x)  q,(t)  +  y  <p;y(x)  qj(t)  1  j 
i  i 

+  [  y  <PiZ(x)  qj(t)  ]  k 

i 


E.  KINETIC  ENERGY 

The  kinetic  energy  of  the  system  can  be  expressed  by  the  summation  of  three 
different  kinetic  energies.  One  is  the  kinetic  energy  of  the  boom  itself,  another  is  the 
kinetic  energy  of  the  tip  mass  and  the  other  is  the  kinetic  energy  of  the  R.F  electronic 
box  attached  to  the  origin  of  the  boom. 
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'  The  kinetic  energy  of  the  boom  itself  is 

t  T  m 

lbm 

1  A* 

—  J0R(x,t)  •  R(x,t)  dm 

(eqn  2.14) 

- 

1  • 

—  J  R(x,t)  •  R(x.t)  pdx 

*0  V 

where 

dm : 

differential  mass  of  the  boom 

dx  : 

differential  length  of  the  boom 

P  : 

mass  per  unit  length 

i  : 

total  length  of  the  boom 

R(x,t) : 

absolute  velocity  of  the  points  on  the  boom 

1  The  kinetic  energy  of  the  tip  mass  is 

l  •  • 

—  M  R(€.t)  •  R(£,t) 

(eqn  2.151 

|  where 

M  : 

tip  mass 

R(£.t>  : 

the  absolute  velocity  at  the  tip  position 

|  The  kinetic  energy  of  the  R.F  electronic  box  is 

t 

i 

1 

Trf=- 

--ir  e2 

2  rzz 

(eqn  2.16) 

where 

»V 

mass  moment  of  inerna  for  the  electronic  box 

time  derivative  of  angular  displacement 

I  Thus  the  total  kinetic  energy  of  the  system  is 

| 

k 
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T  -  T6in  +  Tm  +  Trf 


(eqn  2.17) 


1  .f  •  • 

—  R(x,t)  •  R(x,t)  dm 


+  _L_  M  R(f,t)  •  R(C,t) 


—  ir  e2 

2  rzz 


F.  POTENTIAL  ENERGY 


The  potential  energy  U  of  the  flexible  boom  system  can  be  composed  of  the 
gravitational  potential  energy  U0  due  to  rotation  of  the  system  and  the  strain  energy 

9 

(or  elastic  energy)  L:s  of  the  boom  due  to  deformation.  But  in  our  model  analysis,  we 
exclude  gravitational  acceleration  and  only  consider  the  strain  energy.  The  potential 
energy  was  determined  by  the  work  done  by  the  static  weight  in  the  dellection.  This 
work  is,  of  course  ,  stored  in  the  flexible  flexible  boom  system  as  strain  energy. 

In  this  thesis,  we  apply  mide  summation  method  to  expand  the  deflection  in  terms  of 
the  normal  modes  of  the  system.  The  deflection  of  a  boom  without  any  external  forces, 
satisfies  the  lollowing  equation  of  motion.  [Ref.  7] 


[  ElW'(x.t)  ]"  +  pW(x.t)  =  0 


(eqn  2.  IS) 


where 


mass  per  unit  length 
flexural  rigidity 


and  '  represent  the  derivative  with  respect  to  x.  The  normal  modes  <Pj(x)  of  the  boom 
must  satisfy  the  equation 


EI<p-'(x)  ]"  -  p  <Pj(x)  =  0 


(eqn  2. 19) 


and  its  boundary  conditions.  The  normal  modes  <pj(x)  also  satisfy  the  orthogonality 
relation 


A 

J0<Pj(x)  <Pj(x)  dm  -  0  (for  j  *  i )  (eqn  2.20) 

=  (for  j  =  i ) 

where  M-  is  the  generalized  modal  mass  of  the  i-th  mode. 

As  expressed  in  Appendix  A,  the  deflection  of  the  boom  in  the  equation  A. 4  for 
the  general  form  is 


W(x,t)=  £  <P[(x)  qj(t)  (eqn  2.21) 

i 

and  the  generalized  coordinate  q-(t)  can  be  determined  by  applying  Lagrange's 
equation  after  setting  up  the  kinetic  and  potential  energies. 

Now  the  potential  energy  can  be  expressed  as 


1  r*  7 

u  =  —  J0EIW"2(x,t)  dx 


1  £ 

—I  I  J0  EI<Pi"(x)  <Pf (*)  dx 


i  J 


(eqn  2.22) 


Multiplying  the  both  sides  of  equation  2.19  by  <Pj(x)  and  after  integration  for  the  whole 
boom,  equation  2.19  becomes 


J 


l 

i/Pj(x)  l  EI  <f>i”(x)  1  dx 


(eqn  2.23) 


After  integration  by  parts  and  using  the  boundary  conditions,  equation  2.23  becomes 


'  vwwiwwwwvitfiiniFiuwvTiriif 


<V  J0^x)  <Pj(x)  dm 


From  the  orthogonality  condition,  the  equation  2.20  and  from  equation  2.24,  the  strain 
energy  in  equation  2.22  becomes 


U  *  £  Wj2  Mj  qi{t)2 


(eqn  2.25) 


G.  DERIVATIONS  OF  EQUATIONS  OF  MOTION 

From  the  equation  2.13.  the  dot  product  of  R(x,t)  is 


R(x.f)«R(x,t)  =  [  -  0  £  <prv(x)  q-(t)  +  £  <p-x(x)  qj(t)  J2 


(eqn  2.26) 


+  [  &  Rx(x)+  0  £  «P-X(x>  q.(t)  +  £  <pj-v(x)  q.(t)  l2 


[  £  <PiZ(x)  qs(t>  Y 


®2  (  £  <Piy(x)  q;(t)  ]2  +  [  £  9js(x)  qj(t)  ]2 


-  2  6  £  <p;(x)  q;(t)  £  <p;X(x)  q.(t) 


+  Rx(x)2  02  +  02  (  £  Cp.x(x>  q;( t)  ]2  +■  I  £  (p;-v(x)  q,(t)  l2 

i  l 

4-  2  RX(X)  02  q.(t)  +•  2  0  £  (p.x(x)  q;(t)  £  <P” (x)  q^t) 


2  Rx(x)  0  £  <p;-(\i  *  (t)  +  [  V  (p.z(x)  q.(t)  ]‘ 


=  [  {  £  <P,y(x)  q;(n  ]2  +  {  £  <pjX(x)  q.((t)  }2  4-  {  £  <p;Z(x)  q.(t)  }2 

i  i  i 

+  [  \  £  <P;X( x>  q;(t>  !“  +  (  £  <p.>(x)  q^t)  !2  4-  <  £  (p.zf x )  qit)  !2  ] 


-  2  e  y  fu>‘(x)  qj(t)  V  <p.x(x)  q;(t)  +  Rx2(x)  02 
i  i 

+  2  Rx(x)  e2  X  <P*(x)  q4(t>  +  2  0  y  (p^x)  q4(t)  X  <pjy(x)  qj(t) 
i  i  i 

+  2  Rx(x)  0  X  <Pf(x)  qj(t)  -  02  {  X  <P;Z(x)  q-t(t)  }2 
i  i 

Substituting  equation  2.26  into  the  equation  2.17  and  apply  orthogonal  relationship 
£ 

f  l  <p*(x)  <PjX(x)  +  <p;y(x)  (pj-v(x)  +  (pjz(x)  <p.z(x)  }  dm  (eqn  2.27) 

+  M  [  <p;x(£)  <PjX(C)  +(pj-v(C)  <prv(C)  +  (p.z(C)  tp^C)  }  =  0  (for  i  *  j  ) 

=  Mj  (for  i  =  j  ) 

then  the  kinetic  energy  is 

T  -  — J-  02  X  *,“(*)  M;  +  *4-  X  q5(^2  M.  (eqn  2.2S) 

“  i 

•  _  .  .£ 

-  0  x  X  q;(t>  qj(t)  l  ^  <p;y(x)  <p.x(x>  dm  +  M  <p;>(£)  <p.x(C)  ] 

1  •  i  f  6  ■> 

+  - 0*  f  f  R  “(X)  dm  +■  M  R  ill 

2  -0  x  x 

.£ 

-  v  q-(t)  [  1'  R  »x)  <p.x(x)  dm  ~  MR  (£)  (p.x(E) 

1  «>A  A  1  XI 

i 

•  •  £ 

+-  0  X  X  q;(t)  q.(t )  [  j  <p;x(x)  <pjV(x)  dm  M  (p;xf£)  <p.y(£)  | 
i  ]  O  ' 

•  _  .  .  £ 

+  0  X  q;(t)  (  J  Rx(x)  <piy( x )  dm  +  M  Rx(£)  <p.-v(£)  ] 
i  ^ 

“  02  [  ]'  1  I  <PjZ(x)  qs(t)  }2  +  M  ;  X  <P:'7(C)  q=(t)  l2  ] 


By  applying  the  kinetic  energy  and  potential  energy  expressions  (  eqn  2.25  )  and  (  eqn 
2.28  ),  to  the  Lagrange's  equations  2.7  and  2.8  then  equations  of  motion  will  be 
reduced  as  follows.  The  detailed  derivation  process  can  be  formed  in  Appendix  B. 

•  •  £ 

[  I  q;2(t)  M;  +  f  RX2(X)  dm  +  Vi  Rx2(0  (eqn  2.29) 

i  y 

+  2  I  qj(t)  {  J  £Rx(x)  <PiX(x)  dm  +  M  Rx(0  <p;x(£) } 
i  u 

~  2  (  J  J  I  <PjZ(x)  q}(t)  )2  +  M  (  £  ¥•*(*)  qj(t)  )2  +  Ir  ] 

’  i  i  zz 

•  .  £ 

+  2  0  y  q.(t)  [  qt(t)  Mj  +  f  Rx(x)  <p;x(x)  dm  +  M  Rx(£)  <p;x(C>  1 
i  ^ 

£ 

■  2  r  lj(t)  (  J0<PiZ(x)  <PjZ(x>  dm  -  M  tpjZ( C)  <PjZ(C) }  ] 

£ 

~  y  q*(t)  [  y  q.(t)  {  ^  <px(x)  <p;y( x )  dm  +  VI  (p;X(E)  <pjX(C) 

£ 

”  f0<Piy(x)  <PjX(x)  dm  -  M  (pj'v(C)  <pjX(C)  ) 

.e 

-  JQRX(X)  <p.y(x)dm  -  M  RX(C)  «pi>(E)  J  =  t 

••  •  •  _  c 

<lhW  Mh  +  0  l  Z  <ljW  <  i  <P;X(x)  <Phy(x)  dm  •+■  M  (p.V E)  <phy( E)  (eqn  2.30) 

i 

£ 

“  J0<Piy(x)  <phx(x)  dm  -  VI  <p;y( C)  <phx(£) ) 

£ 

+  JqRx(x)  <phy(x)  dm  +-  Vt  RX(E)  <phy( E)  ] 

•  •  .  £ 

+  29^  qj(t)  [  <PjX(x)  (phy(x)  dm  +  VI  <p;x( E)  <phy(£) 

£ 

-  J(^p.y(x)  (phx(x»  dm  -  VI  (p * y ( E)  <phx(E)  ] 


-  82  [  q,(t)  M,  -  J*  *,(*)**(*>  dm  -  MR,<t)  *,«({) 

+  2  £  q.(t)  {  J  *<pjZ(x)  q>hx(x)  dm  +  M  */(£)  <phz(*)  }] 

+  <°2h  Mh  *  0 

(  h  -  1,  2,  3, . n) 

Now  let's  define  the  following  quantities: 

M0  -  [  £  q.2(f)  M;  +  JqRx2(x)  dm  +  M  Rx2(£) 

+  2  I  q,(t)  {  JlV(x)  <PiX(x)  dm  -t-  M  Rx(£)  ^(f)}  • 

1  * 
t 

~  2  {  f  (  £  <PiZ(x)  q;(t)  )2  +  M  (  V  <p-tz( C)  q.(t)  )2  +■  Ir  ] 
U  i  j  zz 

M0q.  “  2  l  Mi  +  J'jj  R*(x)  <PiX(x>  dm  +  M  Rx(£)  <p,x(«) } 

“  2  X  qj(t)  (  J^«Pjz(x)  <PjZ(x)  dm  -  M  <PjZ(£)  <pjZ(*) }  1 

£ 

Mq  “  “  1 1!  qj(t)  {  j  <P;X(x>  <prvfx>  dm  +  M  <pjx( £)  <PjX(£) 

t  i  U 

e 

-  J0(P;y(x)  <PjX(x)  dm  -  M  <p.y'£)  (PjX(£)  1 

-  J^Rx(x)  <p;y(x)  dm  -  M  Rx(£>  <pjy(f)  1 

i 

•M0h  ■  I Z  q^  i  L<p-,Vx)  <phy(-x)  dm  +  M  <pix(£)  <Phy(fv 

i  ’ 

-  J^<p./(x)  <phx(x)  dm  -  M  <Pjy(£)  <phx(f) } 

£ 

J  .  R  (x)  <phyVx)  dm  +■  M  R  (t>  tphy( £>  ] 


Me<ih  -  2  [  J  *q>.*(x)  q>hy(x)  dm  +  M  q>.x(£)  <f>hy«) 

"  J*<Piy(x)  «phx(x)  dm  -  M  <p;y(C)  <phx(()  ] 

Fch  "  I  <Jj(t)  Mj  -  f  *Rs(x)  <p/(x)  dm  -  MRX(C)  <p.x(C) 

+  2  J  qj(t)  {  J  Jp^x)  <phx(x)  dm  +  M  <p;Z(t)  <phz(C) }] 

^  "  M  %2  %(<> 

Then  equation  2.29  becomes 

M0*8  4-  0  Y  q.(t)  +  Mq  q*(t)  =  t  (eqn  2.31) 

and  quation  2.30  becomes 

M0h  0  +  Me  0  Y  q.(t)  +  Mh  q.(t)  -  F  +  F  *  0 

n  n 


(eqn  2.32) 


III.  COMPUTER  IMPLEMENTATION 


A.  INTRODUCTION 

To  solve  the  equations  of  motion,  two  computer  program.  Dynamic  Simulation 
Language  (  DSL  )  and  NAsa  STRuctural  ANalysis  (  NASTRAN  )  program  are  used. 
NASTRAN  is  a  general  purpose  digital  computer  program  for  the  analysis  of  large 
complex  structures.  [Ref.  S]  This  finite  element  computer  program  was  used  in  the 
modal  analysis  which  determine  the  mode  shapes,  generalized  modal  mass,  generalized 
stiffness,  natural  frequency  of  the  LFMR  model.  Then  these  properties  directly 
inputed  to  the  DSL  program  to  get  a  set  of  solutions. 

DSL  is  an  IBM/VS  FORTRAN-based  simulation  language  for  digital  simulation  of 
continuous  system.  [Ref.  9]  It  is  one  of  the  most  effective  for  the  solution  of 
continuous  modeling  and  simulation  problem  with  computational  power  (automatic 
double  precision  and  accurate  timing),  whether  the  problem  is  time  based  or  not. 

For  the  integration  method  to  solve  simultaneous  nonlinear  second  order  coupled 
ordinary  differential  equations.  Runge-Kutta  method  was  chosen.  Runge-Kutta  fifth 
order  integration  method  (  RK5  )  is  self-starting,  stable  and  automatically  determine 
the  step-size  but  this  needs  excessive  computer  time. 

Two  different  cases  are  analyzed.  One  is  the  LFMR  boom  made  of  Aluminum 
Alloy  and  the  other  is  the  LFMR  boom  made  of  Isotropic  Graphite  Epoxy  composite 
material  (  T300/502S  (  0/90:45/ -45)g  ).  The  mass  distribution  of  the  two  cases  are  the 
same.  Therefore  the  boom  model  of  graphite  epoxy  composite  material  was  suffer  since 
the  mass  density  was  linear. 

B.  MODAL  ANALYSIS 

As  previously  mentioned,  we  consider  our  model  as  a  one  body  system  and  we 
don't  need  any  compatibility  conditions  which  is  necessary  in  multibody  dynamics.  For 
the  modal  analysis  we  equally  divide  the  whole  boom  with  fourteen  grid  points.  Figure 
3.1  shows  Finite  element  model  of  LFMR  boom.  As  the  number  of  grid  points 
increases  we  can  get  more  accurate  results.  But  Degrees  of  Freedom  (  DOF  )  of  the 
system  also  increases.  Fourteen  grid  point  is  sufficient  lor  our  our  analysis. 

In  the  modal  analysis  of  the  boom  .  Modified  Given's  method  (  MGIV  >  was 
applied  and  for  the  purpose  of  simplifying  the  equation  of  motion  wo  normalized  die 


mode  shape  such  that  the  generalized  modal  mass  equal  to  unity.  From  the 
relationship  of  the  generalized  stiffness,  generalized  modal  mass  and  natural  frequency, 
the  generalized  stiffness  Kj  is 

and  reduces  to 

Ki  =  ®j2 

Only  the  first  2  modes  were  needed  in  the  3-dimensional  dynamic  analysis  with 
sufficient  accuracy. 

NASTRAN  program  is  shown  in  Appendix  C  and  some  outputs  are  tabulated  in 
Table  1  and  2  for  3-dimensional  motion.  Figures  3.2  and  3.3  show  the  first  and  the 
second  mode  shapes  for  three  dimensional  motion.  In  Figures  3.2  and  3.3,  the  left  side 
figure  is  the  projection  to  x-z  plane  and  right  side  figure  is  projection  to  v-z  plane. 

The  first  mode  shape  shows  there  is  no  in-plane  vibration  and  the  second  mode 
shape  shows  there  is  no  out-of-plane  vibration. 

TABLE  1 

REAL  EIGENVALUES  OF  ALUMINUM  ALLOY  (3D) 


mode 

no. 

radians 

cvcles 

generalized 

mass(M;) 

generalized 
stillness(  Kd 

1 

3.3845 15E  + 00 

5.386623E  — 01 

1.000000E  +  00 

1.145494E  +  01 

2 

3.568821  E  + 00 

5.679954E-01 

— 

1 .000000 E  -t-  00 

1.273648E  +  01 

C.  DSL  PROGRAMING 

DSL  was  implemented  for  the  solution  of  a  set  of  simultaneous,  nonlinear, 
second  order,  ordinary  differential  equations.  DSL  oilers  a  new  simulation  tool  that 
speeds  problem  analysis  and  solution  for  a  wide  variety  of  application.  It  is  adopted  for 
simulation  because  it  requires  less  skillful  time  for  problem  analysis  and  provides 
quicker,  more  comprehensive  plotting  through  G RAFAEL  and  sophiscated  line  or 
print  plot.  These  advantages  can  translate  directly  into  higher  productivity  and  cost 
savings. 
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TABLE  2 

REAL  EIGENVALUES  OF  COMPOSITE  MATERIAL  (3D) 


mode 

no. 

radians 

(O- 

i 

cycles 

0)' 

i 

generalized 

stiffnessfKj) 

1 

4.287485E  +  00 

6.823744E-01 

1.000000E  +  00 

1.S38253E  -‘-•*1 

2 

4.461334E  +  00 

7.100433E-01 

1.000000E  +  00 

1.990350E-01 

DSL  is  a  high-level  continuous  simulation  language  which  incorporates  VS 
FORTRAN  as  a  subset.  Because  of  its  tremendous  flexibility,  DSL  facilitates  the 
solution  of  nearly  all  prlblems  involving  time-dependent  differential  equations.  Thus 
DSL  readily  assists  in  the  dynamic  analysis  of  transient  behavior  of  dynamic  systems. 
Also  DSL  is  easily  learned  and  applied  to  many  problems  in  science,  engineering, 
mathematics  and  management.  Coding  is  simple,  execution  is  rapid  and  results  can  be 
displayed  graphically.  For  any  one  involved  in  simulation  modeling  DSL  otfers 
increased  power  for  faster  problem  solving  and  the  user  choice  of  nine  integration 
method  (  fixed-step,  variable-step,  variable-step  &  variable-order  method  ). 

In  our  model  analysis.  Runge-Kutta  fifth  order  method  (  RK5  >  was  used  because 
it  is  self-starting,  stable  and  provides  good  accuracy.  To  code  the  equations  of  motion, 
the  equation  2.31  and  2.32  have  to  be  rewritten  as  follows: 


VIq  0  +  Mq  q.(t)  =  r  —  Mq  0  V  q.(t)  (eun  3.1) 

and  equation  2.30  becomes 

M0h  0 +  Mh  =  Fc.  “  Feh  ~  0  I  ^ 

n  n 

(  h  -  l.  2.  3 . n  ) 

Two  ways  solving  these  systems  of  equations  are  shown.  One  uses  matrix  algebra  and 
utilizes  subroutines  from  the  library  LINPACK,  the  other  uses  subroutine  REGULRA 
provided  by  DSL.  For  our  system  analysis  we  select  matrix  algebra  because  it  is  easy 
to  use  for  iarge  numbers  of  variables. 


Using  matrix  notation,  the  equations  3.i  and  3.2  may  be  rewritten  as 

I M  1  {  X  >  -  {  f } 


where  the  matrices  M,  X,  and  f  are  defined  as  follows: 


and  by  successive  integration  of  {  X  },  the  solution  vector  may  be  obtained.  Since  the 

•  • 

elements  of  the  [  M  ]  and  {  f  }  matrices  are  time  dependent,  the  {  X  }  vector  must  be 
computed  at  each  integration  step.  A  matrix  decomposition  subroutine  (  DGEFA  ) 
using  Gaussian  elimination  and  a  subroutine  which  uses  this  decomposition  to  solve  a 

matrix  equation  (  DGESL  )  are  called  directly  from  the  LINPACK  library.  Because 

•  • 

DGESL  returns  the  solution  vector  {  X  }  in  the  right-hand-side  vector  f,  the  vector 

•  • 

name  {  X  }  does  not  appear  explicitly  in  the  program. 

The  coefficients  of  [  M  ]  and  (  f }  matrices  are  calculated  in  each  time  step.  9,  qj 

(  i  =  1.  2,  3, . .  n  ),  deflection  and  siope  at  the  tip  position  in  each  direction  were 

calculated  at  each  time  interval  of  interest.  The  graphs  of  these  variables  vs.  time  were 
also  obtained.  The  computer  program  for  dynamic  analysis  is  coded  in  Appendix  D. 


IV.  RESULTS  AND  DISCUSSIONS 


Computer  simulation  was  done  in  four  areas  for  double  link  boom  system  in 
three  dimensional  motion  to  investigate  the  equilibrium  configuration  and  the  vibration 
amplitude  of  the  tip  position.  The  effects  of  the  spin  up  procedure  on  the  pointing 
error  of  the  reflector  is  studied;  the  first  analysis  was  the  comparison  of  the  three 
different  torque  histories  by  maintaining  magnitude  of  torque  until  rotating  speed  of 
boom  reaches  to  15  r.p.m.  Secondly  the  magnitude  of  appllied  torque  was  changed  in 
three  cases  for  one  of  three  torque  applying  methods  above.  The  effects  of  structural 
damping  of  the  boom  on  the  settling  time  is  also  studied  by  changing  the  modal 
damping  values.  The  magnitude  of  the  deflection  and  slope  at  the  equilibrium 
configuration  at  three  different  rotating  speeds  are  investigated  by  simulating  free 
motions  of  the  system  with  the  three  different  initial  rotating  speed  and  undeformed 
configuration.  Also  some  comparisons  were  made  for  double  link,  flexible  boom  system 
in  planar  motion  in  Appendix  B. 

A.  MATERIAL  PROPERTIES  AND  PARAMETERS 

For  the  comparison  of  the  computer  simulation  results,  two  kinds  of  material  were 
chosen.  One  is  Aluminum  alloy  6061-T6  (99%  Al-1%  Mg)  and  the  other  is  an 
Isotropic  Graphite-Epoxy  Composite  material  (T300'5208  (0:90/45  -45)g).  [Ref.  10] 
Same  mass,  same  length  of  each  link,  same  outer  diameter  was  used  for  these  two 
materials.  Size  of  the  electronic  box  is  1  ft.  cube.  The  properties  of  the  materials  and 
some  geometric  parameters  are  given  in  Table  3  and  Table  4  respectivly. 

B.  EFFECTS  OF  TORQUE  APPLYING  PROCEDURE 

Three  cases  of  torque  applying  method  were  chosen  without  damping  for 
comparison.  To  maintain  15  r.p.m  of  rotating  speed,  w’e  select  the  maximum  magnitude 
of  torque  is  10  Ibs-in.  This  value  may  be  appropriate  to  maintain  15  r.p.m  in  20  to  30 
seconds. 

In  the  first  case,  10  lbs-in  torque  was  abruptly  applied  from  the  begimng.  This 
method  is  unlikely  in  the  real  situation  but  we  chose  this  case  for  comparison  purpose. 
In  the  second  case,  torque  was  parabolically  applied  until  reaching  15  r.p.m  then  cut 
off  abruptly  to  0.  Finally  torque  was  applied  same  manner  as  case  two  but  linearly 
decreased  to  0.  Figure  4.1  shows  three  different  spin-up  schemes. 
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TABLE  3 

BOOM  PROPERTIES 


PROPERTY 

ALUMINUM  ALLOY 

COMPOSITE 

specific  weight  (  y  ) 

0.09S0  lb  in3 

0.058S  lb  in3 

modulus  of  elasticity  (  E  ) 

10.  IE  +  06  psi 

10. IE +  06  psi 

modulus  of  rigidity  (  G  ) 

3.7E  +  06  psi 

4.1E  +  06  psi 

poisson's  ratio  (  v  ) 

0.360 

0.250 

outer  diameter  (  r  j  ) 

3.0  in 

3.0  in 

inner  diameter  (  r->  ) 

2.7312  in 

2.5362  in 

thickness (  t  ) 

0.1344  in 

0.2319  in 

cross  area  (  A  ) 

1 .2 1 0 1 2 

2.016S  in2 

mass  per  unit  length  (  p  ) 

3. 0690 E -04  lb  in 

3.0690E-04  lb.  in 

area  moment  of  inertia  (  I  ) 

1.2447  in4 

1.9451  in4 

polar  moment  of  inertia  (  J  ) 

2. 4894  in4 

3.S902  in4 

In  these  three  cases  angular  displacement  8  was  changed  parabolically  as  shown 
in  Figure  4.2  while  torque  is  applied,  then  linearly  increased  as  torque  was  cut  olF.  The 
cases  3  and  C  are  almost  indentical  and  overlapped  in  the  Figure  4.2.  Angular 

displacement  was  equally  increased  with  time  lor  aluminum  ailoy  and  composite 

• 

material  because  total  system  mass  is  same.  As  shown  in  Figure  4.3,  angular  velocity  8 
varied  linearly  with  some  oscillation  in  the  begimng  for  the  first  case  because  torque 
was  suddenly  applied,  then  it  gradually  decreased  until  reaching  15  r.p.m.  Alter  cutting 
olF  applied  torque,  8  oscillates  periodically  because  the  boom  deflects  to  the  positive 
direction  and  the  concentrated  mass  center  moves  to  the  rotating  axe',  this  means 
radius  of  rotation  decreased  and  moment  of  inertia  also  decreased.  From  the 
conservation  of  angular  momentum. 

Hq  =  m  r2  8  =  constant 
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TABLE  4 

GEOMETRIC  PARAMETERS 


PARAMETER 

— 

VALUE 

Cj  ;  length  of  lower  boom 

168  in 

Cs  ;  length  of  upper  boom 

144  in 

a  ;  angle  between  Cj  and  x-axes 

-1 

O 

O 

P  ;  angle  between  f|  and  f-> 

126° 

M  ;  tip  mass 

37.5  lbs 

rr^  ;  mass  of  R.F  electronic  box 

•  50  lbs 

*  *> 

So  as  the  tip  position  changes,  0  varies  with  reciprocal  of  r“.  Therefore  angular 
velocity  oscillates  periodically. 

In  case  two  and  three,  0  increased  smoothly  without  oscillation  until  reaching  15 
r.p.m  then  it  had  relatively  small  oscillation  as  soon  as  torque  was  removed  gradually. 
Figure  4.5  shows  the  angular  velocity  of  the  Aluminum  Alloy  boom  in  magnified  scale. 
For  the  composite  material  as  shown  in  Figure  4.4,  magnitude  of  oscillation  was  much 
smaller  than  that  of  aluminum  alloy  so  it  looks  no  oscillation.  But  magnified  angular 
velocity  of  constant  part  shows  obvious  oscillation  in  Figure  4.6.  1  ires  4."  -  4.10 
show  the  variations  of  the  generalized  coordinates  qj  and  q,  respectively. 

During  rotation  dellection  in  x-direction  was  dominated  and  Jetlection  center 
readies  to  its  equilibrium  position  in  x,  /-direction  then  oscillates  harmonically  as 
shown  in  Figures  4.11  and  4.12  but  displacement  in  y-direction  was  rapidly  increased  at 
the  begining  then  decreased  gradually  to  0  and  as  soon  as  applied  torque  was  removed 
it  oscillates  harmonically  as  shown  m  Figures  4.13  and  Figure  4.14. 

Characteristics  of  displacement  and  slope  at  the  tip  position  is  very  similar  to  its 
generalized  coordinates,  and  the  generalized  coordinates  q,  in  Figures  4.4  and  4.F'  was 
dominated  for  the  system  dellection.  Displacement  in  /-direction  was  shown  in  Figure 
4.15  and  16  for  two  material  booms.  Displacement  and  slope  change  as  shown  m 
Figures  4.17  -  4.22  was  apparently  dilferent  from  case  A  and  case  B  out  alter  cutting 
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off  the  torque  its  characteristics  of  deformation  was  not  much  different.  For  the  case 
three  deflection  and  slope  increases  smoothly  while  torque  is  applying  then  as  torque 
was  removed  it  reaches  its  equilibrium  position  and  oscillates  with  small  fluctuation 
error. 

Figures  4.11  -4.24  show’  displacement  and  slope  at  the  tip  position  in  each 
direction  for  aluminum  alloy  boom  and  composite  material  boom. 

As  we  mentioned  previously,  we  are  mainly  interested  in  the  pointing  error  at  the 
tip  position  in  both  elevation  and  azimuth.  For  the  composite  material  boom,  all  3 
cases  were  satisfied  the  pointing  error  requirement  (  elevation  and  azimuth  angle  error  : 
0.0328°  ).  But  for  the  Aluminum  Alloy  boom,  torque  applying  case  A  and  case  B  were 
not  satisfied  the  elevation  angle  error  requirement.  Figures  4.19-4.22  show  the 
characteristics  of  elevation  and  azimuth  angle  variation,  and  in  Figures  4.23  and  4.24, 
the  effects  of  torque  applying  procedure  are  compared.  From  these  results,  pointing 
error  at  the  equilibrium  position  depends  on  torque  removing  procedure.  As  we  see  in 
Figures  4.23  and  4.24.  pointing  error  in  case  A  and  case  B  is  almost  same  but  it  is  very 
different  in  case  C  and  pointing  error  varies  more  sensitively  with  flexible  material  as 
torque  removing  procedure  varies.  Consequently,  if  we  apply  the  torque  more  gently 
and  remove  slowly  wath  sufficient  time,  the  pointing  error  both  deflection  and  the  slope 
can  be  reduced  to  much  smaller  values. 

C.  EFFECTS  OF  CHANGING  MAGNITUDE  OF  TORQUE 

In  this  section  we  investigate  the  effects  of  the  maximum  magnitude  of  applied 
torque.  For  the  simple  comparison,  we  select  second  case  of  torque  applying  method 
used  in  previous  section.  Magnitudes  of  torque  chosen  are  10  lbs-in.  20  Ibs-in  and  40 
lbs-in  as  shown  in  Figure  4.25.  As  the  magnitude  of  applied  torque  increases  4  times, 
applying  time  was  reduced  to  almost  one  third.  All  characteristics  are  same  as  the 
second  case  of  previous  section  for  3  cases  shown  in  Figures  4.26-4.36  but  magnitude 
of  pointing  error  increased  linearly  as  the  maximum  torque  increases  as  shown  in 
Figures  4.37  and  4.38.  Figures  4.26-4.43  show'  angular  displacement,  angular  velocity, 
generalized  coordinates,  displacement  and  pointing  error  for  the  boom  made  of 
Aluminum  Alloy. 

D.  EFFECTS  OF  DAMPING  COEFFICIENT 

In  this  section,  we  will  investigate  the  damping  effects  for  settling  down  the 
vibration  of  the  boom.  As  we  mentioned  before,  in  the  3  cases  of  torque  applying 


procedure,  the  second  case  did  not  satisfy  the  pointing  error  requirement  in  elevation 
angle.  We  will  not  mention  about  the  first  case  because  it  is  unlikely  to  be  used  in  real 
situations.  For  the  comparison  purpose,  we  arbitrary  choose  3  modal  damping 
coefficient  such  as  0.2  %,  0.5  %  and  1.0  %.  To  see  the  results  more  clearly,  we 
simulate  the  model  in  200  seconds.  Figure  4.39  shows  the  spin- up  procedure  used  in 
this  analysis.  Figure  4.40  shows  vibration  settling  down  ratio  during  the  first  200 
seconds.  In  this  Figure,  it  is  hard  to  recognize  the  effects  of  all  3  cases  but  in  Figure 
4.41,  we  can  clearly  see  the  vibration  settling  down  for  each  case. 

Two  trials  were  made  for  investigation.  Firstly,  we  tried  to  find  how  much 
vibration  was  settled  down  in  200  seconds  for  each  damping  coefficient.  Secondly,  we 
found  how  much  time  was  required  to  meet  the  pointing  error  requirement  for  each 
case.  Figure  4.42  shows  elevation  angle  pointing  error  decreasing  ratio  in  200  seconds. 
Each  value  stands  for  magnitude  of  pointing  error  ratio  to  its  initial  value  which  is  the 
magnitude  of  pointing  error  without  damping.  As  shown  in  Figure  4.42  pointing  error- 
decreased  exponentially  as  damping  coefficient  increases  linearly.  Figure  4.43  shows 
desired  time  to  meet  pointing  error  requirement  for  each  damping  coefficient.  As 
damping  coefficient  increases  linearly,  desired  time  decreased  exponentially. 

E.  THE  EQUILIBRIUM  CONFIGURATION  AT  DIFFERENT  ROTATING 

SPEEDS 

Initial  conditions  were  given  for  3  cases  such  as  5  r.p.m.  10  r.p.m  and  15  r.p.m 
with  damping  coefficient  of  2  %.  From  this  result,  we  could  investigate  the  magnitude 
of  deflection  in  each  direction  and  slope  at  the  tip  position.  Figure  4.45  shows  angular 
velocity  (  rotating  speed  )  variation  with  time.  As  the  initial  rotating  speed  increases, 
magnitude  of  oscillation  was  aiso  increased  and  the  rotating  speed  at  the  equilibrium 
state  also  increased.  As  we  mentioned  in  section  B.  angular  momentum  is  constant  in 
the  system  since  no  external  torque  is  applying.  Therefore  as  the  boom  defects  to  the 
positive  x-direction,  radius  of  rotation  decreases  so  angular  velocity  increases 
consequently.  Figure  4.45  shows  elevation  angle  change  and  its  magnitude  of 
fluctuation  change.  As  the  magnitude  of  initial  speed  increases  linearly,  elevation  angie 
change  and  magnitude  of  fluctuation  increases  exponentially.  This  result  is  shown  in 
Figure  4.47.  Magnitude  of  deflection  in  x-direction  and  z-direction  was  also  increased 
as  the  rotating  speed  increased  as  shown  in  Figures  4.4”  and  4.48  respectively.  This 
effects  was  more  sensitive  for  the  flexible  material  boom.  The  oscillation  of  azimuth 
angle  centered  at  •». 
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FIGURES 


Figure  4.2  Angular  displacement  vs.  time 
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Figure  4.3  Angular  velocity  vs.  time  t  AL.  ). 


Figure  4.6  Magnified  angular  velocity  vs.  time  (  COM.  ). 
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Figure  4.9  Second  mode  generalized  displacement  vs.  time  <  AL.  ). 
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Figure  4.12  Displacement  in  x-uircction  vs.  time  t  COM.  ). 
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Figure  4.14  Displacement  in  y-direction  vs.  time  (  COM.  ). 
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ELEVATION  ANGLE  ERROR  RATIO 


Figure  4.23  Pointing  error  in  elevation  angle  vs.  torque  applying  procedure. 
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Figure  4.26  Angular  displacement  vs.  time  (  AL.  ). 
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Figure  4.37  Pointing  error  change  in  elevation  angle  vs.  magnitude  of  torque  \  AL.  i. 
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Figure  4.39  Applied  torque  vs.  time  with  damping  (  AL.  ). 


ERROR  DECREASING  RATIO  IN  200 


TIME  NEEDED  FOR  REQUIRMENT 


Figure  4.43  Desired  time  to  meet  elevation  angle  error  requirement  vs.  damping  coeificient. 
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Figure  4.44  Initial  angular  velocities  vs.  time  with  damping  (  AL.  ). 
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V.  CONCLUSIONS  AND  RECOMMENDATIONS 


A.  CONCLUSIONS 

The  purposes  of  this  research  are  to  investigate  the  the  effects  of  the  flexibility  of 
the  LFMR  boom  on  pointing  error  of  the  reflector  in  elevation  and  azimuth  angie  and 
to  identify  the  parameters  which  are  important  in  pointing  error.  In  this  study,  the 
effectiveness  of  the  application  of  Lagrange's  method  and  mode  superposition 
techniques,  computer  simulation  techniques  are  also  investigated. 

The  results  indicated  that  the  Lagrange's  method  and  mode  superposition 
technique  was  very  effective  for  this  double  link  boom  model  analysis.  The  effects  of 
flexibility  of  boom  was  sufficiently  expressed  with  few  mode  and  the  boom  flexibility 
was  very  affective  in  pointing  error.  As  the  material  becomes  stiffer  with  light  weight, 
its  pointing  error  becomes  smaller.  The  torque  applying  and  removing  procedure  is 
very  important  to  the  magnitude  of  pointing  error.  If  we  apply  the  torque  gradually 
until  reaching  desired  rotating  speed  and  remove  slowly  with  sufficient  time,  we  can 
reduce  the  pointing  error  to  much  smaller  values  within  requirement.  The  vibration 
settling  time  decreases  exponentially  as  modal  damping  coefficient  increases  and  the 
pointing  error  is  linearly  dependent  on  the  magnitude  of  applied  torque.  The  deflection 
in  each  direction  and  elevation  angle  change  in  equilibrium  condition  increase 
exponentially  as  the  rotating  speed  increases.  Sensitivity  of  deflection  depends  on 
flexibility  of  the  boom. 

B.  RECOMMENDATIONS 

In  this  research,  we  regarded  the  deployable  reflector  as  a  concentrated  tip  mass. 
In  this  case,  we  don  t  need  consideration  of  the  flexibility  of  the  reflector  but  in  actual 
model  the  flexibility  of  the  reflector  will  affect  the  accuracy  of  the  pointing  error, 
therefore  flexibility  of  the  deployable  reflector  has  to  be  considered  for  the  future  work. 

As  we  have  seen  in  the  comparison  of  two  material  boom,  the  flexibility  of  boom 
was  very  affective  in  the  pointing  error.  Stiller  material  with  light  weight  will  play 
important  role  in  the  reduction  of  pointing  error  of  the  flexible  boom. 


APPENDIX  A 

DERIVATIONS  OF  THE  EQUATIONS  OF  MOTION 

In  this  Appendix,  we  develop  simplified  dynamic  equations  of  motion  for  2  cases 
of  the  boom  system  with  tip  mass. 

The  large  motion  caused  by  rotation  and  the  small  motion  created  by  elastic 
deformation  will  be  expressed  by  generalized  coordinates  then  apply  Lagrange's 
equations  to  develop  the  equations  of  motion.  In  these  analysis,  the  local  rotary 
inertia  and  shear  deformation  of  booms  are  neglected. 

From  the  equations  2.7  and  2.S  Lagrange's  equations  are 


d  dT  dT  dL 

—  r  ___  i - +  - 

dt  CV  <30  C0 


(eqn  A.l) 


d  <3T  ^  _  aT  t  au 
dt 

(  h  =  l,  2,  3 . n  ) 


(eqn  A.2) 


where 

0  :  the  generalized  coordinates  of  the  large  motion 

qh  :  the  generalized  coordinates  of  the  small  motion 

t :  applied  torque  to  the  system 

n  :  number  of  modes  (or  number  of  degrees  of  freedom) 


1.  SINGLE  LINK  BOOM  SYSTEM  IN  PLANAR  MOTION 
a.  Geometry  of  the  system 

In  model  analysis,  since  the  extension  deformation  is  negligible,  only  the 
bending  deformation  is  considered  in  the  analysis. 


Figure  A.l  Parameters  of  the  single  link  boom  system. 


M  : 

tip  mass 

"V 

mass  of  electronics  box 

t: 

length  of  the  link 

t  : 

applied  torque 

6(t): 

• 

angular  displacement 

0(0  : 

angular  velocity 

R,(x)  : 

position  vector  in  local  x-direction 

W(x,t)  : 

deformation  vector  of  boom 

R(x.t)  : 

position  vector  of  the  point  on  the  boom  after  deformation 

i  : 

unit  vector  of  local  \-Jircction 

j : 

unit  vector  of  local  y-direction 
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k 


ftt 


unit  vector  of  local  z-direction 
unit  vector  of  global  X-direction 
j0  :  unit  vector  of  global  Y-direction 

kf, :  unit  vector  of  global  Z-direction 


b.  Position  and  velocity 

During  rotation  the  boom  deforms  and  the  positions  of  any  point  in  the 
system  can  be  expressed  the  vector  sum  of  a  vector  Rx(x)  from  the  origin  to  x  and  a 
vector  W(x,t)  caused  by  deformation.  Then  the  position  vector  R(x,t)  is  expressed  as 

R(x,t)  =  Rx(x)  +  W(x.t)  (eqn  A. 3) 

where  Rx(x)  is  only  x  dependent  variable 
so 


Rx(x.)  =  Rx(x)  i 

and  VV(x.t)  is  deformation  obtained  from  modal  summation  method 
then 


W(x.t)  »  V  <Pj(x)  q;(t) 
i 


<eqn  2.21 ) 


where 

(Pj(x)  is  i-th  mode  shape  function 

qt( t )  is  i-th  mode  generalized  coordinates 


If  we  consider  the  Jcformation  of  the  boom  consists  of  translation  and  extension,  then 
equation  2.21  becomes 

W(\,t)  =  ^  [  (pT'x)  qdt)  i  +■  <p.v(xt  qa.t)  j  j  <eqn  A.-ri 

i 

where 
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vvN.-  -'■V-.V. 

VA  A  A.'.  aAx.*  m.  V  <  1 . 


q>jX(x)  i  is  i-ih  mode  shape  function  in  extension 
9jy(x)  j  is  i-th  mode  shape  function  in  translation 

From  equation  2.21  and  A.4  the  position  of  the  point  on  the  boom  is 


i  u  M>.i '.Mil 


R(x,t)  *  Rx(x)  i  +  V  [  <pjX(x)  i  +  <p,V(x)  j  ]  q.(t) 


(eqn  A. 5) 


Now  the  velocities  of  the  point  on  the  boom  is  necessary  to  formulate  the 
kinetic  energy  to  apply  Lagrange's  equation.  The  velocity  is  obtained  by  simply 
differentiating  the  equation  A .5 


R(x,t)  =  Rx(x)  i  +  y  l  <p.(x)  i  +  q>j-v(x)  j  ]  q.(t) 


+  V  [  <p*(x)  i  +  <p.-v(x)  j  ]  qj(t') 


*  [Rx  +  X  <J>iX(x)  qs(t>  ]  i  +  1  <Piyfx)  q.(t)  j 
i  i 

+  1  [  <PiX(x)  qj(t)  i  +<piy(x»qi(t)  j  j 


But  time  derivative  of  unit  vector  i  and  j  are 

i  =  0  k  *  i  =  9  j 
j  =  0k  x  j  =  -  0  i 
k  =  0  k  x  k  =  0 

substitute  these  quantities  to  the  equation  A. 6  and  simplify 
then 

•  •  • 

R(.x.t)  *  {  Rx  *  ^  <IVx(x>  qjtt )  j  0  j  -  O  X  Vj-v(x)  qjtt)  i 
i  i 

■+•  ^  (pAlx)  q.(t)  i  -  I<Pi>(x )  t )  j 


(  X  <PjX(x)  qs<t>  ~  0  1  <pj-v(x)  qjin  ]  i 


(eqn  A. 6) 


i  eqn  A.") 


-  -  o 


+  l  Rx(M  5  +  ©  T  <PiX(x>  q^t)  +  y  if.y(x)  qj(t)  ]  j 
i  i 

c.  kinetic  energy  and  potential  energy 

Recall  the  equation  A.7  and  drop  the  all  terms  concerning  tp^x)  and  express 
<p;?(x)  simply  q>.(x) 
then  equation  A.7  reduces  to 


R(x.t)  -  Rx  0  j  -  0  V  <pj(x)  q j ( t )  i  +-  1  <p;y(x)  q^t)  ]  j 


(eqn  A.S) 


=  -  6  1  <pi(x)qi(t>  i  +  [  Rs(x)  0*1  <Pk(x>  q.(t)  ]  j 


From  the  equation  2.17,  the  kinetic  energy  of  the  whole  system  is 


T  “  Tbm  +  +  Trf 


1  A  •  • 

=  - |  R(x.t)  •  R(x,t)  dm 

">  Jn 


+  —  M  R(f.t)  •  R(£,t) 


(eqn  A. 9) 


•>  r 


L  0- 


•  •  •  ^  ^ 

R(x,t>«R(x,t)  =  [  —  0  1  x )  q.(t)  ]"  +•  1  Rr(x)  0  +  1  <p(x)  qj( t >  )“  ieqn  A  1'" 

i  i 

=  02  [  1  ^(x)  Mji.t)  ]2  +  i  Rx( x )  0)  2 

i 

+  2  RX(X)  0  1  <p((X)  qk( t >  -  [  1  <pj(X»  q.i  t )  i2 
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UNCLASSIFIED 


F/0  22/1 


1 Z  9/x)  qj(t)  l2  -  X  <Pi2(*)  q2(0  +  2  V  <Pj(x)  <pj(x)  q^tiq-ft) 
i  i  i’j 

l  Z  <Pj(x)  qj(t)  ]2  -  Z  <Pi2(x)  q2(0  +  2  £  9j(x)  <pj(x)  q,(t)q.(t) 
i  i  i’j 

(for  i  *  j  *  0  ) 

Substituting  into  equation  A.  10 
then 

R(x,t)«R(x.t)  =  82  [  y  <pf(x)  q;2(t)  +  2  V  ^(x)  *.(x)  q.(t)  q.(t)  (eqn  A.l  1) 

i  i’j 

+  Rx2(x)  92  ]  +  2  Rx{x)  0  ^  <p;(x)  qj(t) 

i 

+  Z  <Pi2(x)  v(t)  +  2  Z  <Pi(x>  <pj(x)  q-.Ct)  q/t) 

i  i’j 

(for  i  *  j  ) 

Now  equation  A.9  can  be  rewritten  as 


02  Z  q2(t)  (  <Pj2(x>  dm  +  M  tpftt)  ] 


C 

+  0“  Z  q^1)  q^)  (  J  <P,(x>  <Pj(x)  dm  +  M  q>.(f)  *.(t)  ] 
i’j  u 


+  -J-  Z  v(t)  i 1,Si2fx)  dm  +  M  <p-,2fe)  i 

-  I 


+  Z  q,(l)  qj(U  (  J  <Pi(x>  (p^x)  dm  +  M  <pj(C)  9j(C)  1 

l-j  u 


+  -y-  02  (  RX2(X)  dm  -  M  RA2(£)  j 


£ 


+  0  Z  qs<t)  [  JQRx(x)  <Pi(x)  dm  +  M  RX(C)  tpj(t)  ] 
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(eqn  A.  1 2) 


VuV 


(for  i  *  j  ) 


+  -i-Ir  82 

2  rzz 


From  orthogonality  relationship 


.£ 


JQ*Pj(x)  <pj(x)  dm  +  M  <p;(x)  <pj(£)  *  0  (for  i  *  j  ) 


(eqn  A.  13) 


=  Mj  (for  i  *  j  ) 


where  M.  is  the  i-th  mode  generalized  mass,  then  the  kinetic  energy  of  the  system  can 
be  simplified  as  follows 


1  ^  ,  i  ^  l  ,  i  •  ^ 

—  ^  1  V(t)  M.  +  —  X  qf(t)  M; 

*•  l  “  l 


(eqn  A.  14) 


—  Ir  02 

2  rzz 


—  02  [  /  R  (x)  dm  +  \1  R  2(£) } 
/  *  * 


£ 


+  0  I  qs(t)  [  JQRx(x)  <P(x)  dm  +  M  Rx(£>  *<*)  1 


From  the  equation  2.25,  potential  energy  is 


U  *  — —  y  to2  M.  q.2(t) 

2  -r  1  i  n  x 


(eqn  A.  15) 


d.  Lagrange's  equations 

Substitute  equations  A.  14  and  A.  15  into  equations  2.19  and  2.20  to  apply 
Lagangie's  equation, 
then 


JT 


cO 


=  0 
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-gg-  ~  » I  V(t)  M,  +  8  [  ^  Rx2(x)  dm  +  M  Rx2«)  ] 


+  I  8|(‘>  1  Sg  *,«  *i<*)  <im  +  M  Rx(t)  9,(1)  ) 


+  I.  8 

*zz 


e  ( I  qr(t)  Mj  +  J0R/(x)  dm  +  M  Rx(t)  +  1^  ] 

t 

+  Z  ^  1  J0  Rx<X)  *,(*)  dm  +  M  Rx(£)  <Pi(£)  ] 


"37 [ 4r 1 '  ® 1  ?  qi:<t)  M1  +  ^oR'2<x)  dm  +  M  R»2<e)  +  lh 

4-2  0  5;  q.(t)  q.(t)  M; 


Z  qj(t)  ( 10  Rx(x)  %(K)  dm  4-  M  RX(C)  <p;(C)  ] 


dV 

6U  d 

<30 

<3(5  dt 

dT 

d% 

’  %(tl  M„ 

dT 

• 

• 

*% 

=  qh(t)  Mh  +  0 

d 

61 

— 1 1 

iz 1  =  qh(t)  M> 

1  -  0 


..  e 

01  v 


REPRODUCED  AT  GOVERNMENT  EXPENSE, 


-  Mh  %w 


dL  d  dU 

-TT-  =  — [  —  ]  -  0 
5qh  lit  cqh 


From  all  these  quantities,  equation  A.l  becomes 


*0  [  V  q  2{t)  \|  +  \  R 2(X)dm  +  M  R  2(C)  +  I.  ] 

~  1  1  *0  *  *  lzz 


(eqn  A.  16) 


+  2  0  v  q.(u  q  ( t )  M.  +  v  q.(t)  [  f  R  (\)  <p  ( x )  dm  +  M  R  (£)  <p(t)  ]  =  T 

11  I  *1  »  A  A  1  A 


and  equation  A. 2  becomes 


0  I  LRx(x)  dm  +  M  Rx(£>  Vh(*)  1  +  V1)  Mh 


(eqn  A.  17) 


-  e2  qh(t)  Mh  +  toh2  Mh  qh(t)  =  ~  2  U  (0h  Mh  qh(t) 


(  h  -  1.  2.  3. 


In  the  equation  A.  16  and  A.  17,  if  we  select  n  numbers  degree  of  freedom  of 
the  system,  that  is,  if  we  select  n  numbers  of  mode  shape,  we  could  have  n+  1  numbers 
of  equations. 


. V.v  V'.v.* 

■  V  ■.  K  M.  JL  M.  «. 


REPRODUCED  AT  GOVERNMENT  EXPENSE,’ 


DOUBLE  LINK  BOOM  SYSTEM  IN  PLANAR  MOTION 
a.  Geometry  of  the  system 


Figure  A. 2  Parameters  of  the  double  link  boom  system  in  planar  motion. 


R0(x) 


tip  mass 

mass  of  electronics  box 
length  of  lower  boom 
length  of  upper  boom 
angle  between  two  links 
applied  torque 
angular  displacement 
angular  velocity- 

position  vector  of  the  point  on  the  boom  in  the  local  x-dircction 


REPRODUCED  AT  GOVERNMENT  EXPENSE,' 


R(x,t) : 

position  vector  of  the  point  on  the  boom  after  deformation 

W(x.t)  : 

deformation  vector  of  boom 

i  : 

unit  vector  of  local  x-direction 

j  = 

unit  vector  of  local  y-direction 

k  : 

unit  vector  of  local  z-direction 

'o  : 

unit  vector  of  global  X-direction 

io : 

unit  vector  of  global  ^'-direction 

kg  : 

unit  vector  of  global  Z-direction 

b.  position  and  velocity 

In  this  double  link  system,  we  have  to  consider  the  extension  deformation  as 
well  as  the  bending  deformation. 

The  position  vector  of  a  point  on  the  boom  was  expressed  as  a  sum  of  Rx(x) 
and  W(x.t). 

Then 

(eqn  2.9) 


(eqn  A. IS) 

% 

From  equation  A.4  and  A.  18  the  position  vector  of  the  point  on  the  boom  is 

R(x,t)  =  Rx(x)  i  +  Ry(x)  j  +  X  l  tPj(x)  i+  q»jy(x)  j  ]  q^t)  (eqn  A.19) 

i 

Now  by  differentiating  the  equation  A.19,  we  obtain  the  velocity  of  the  point  on  the 
boom. 

R(x.t)  =  Rx(x)  i  +  Ry(x)  j  +  V  |  <p.(X)  j  +  <p.>(x)  j  )  qj(t)  (eqn  A.2U) 

i 

+  X  (  <PjX(x)  i  +  <pjy(x>  j  1  q4(t) 


R(x,t)  =  R0(x)  +  W(x,t) 
where 

R0(x)  =  Rx(x)  +  Ry(x) 

=  Rx(x)  i  +  Ry(x)  j 
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REPRODUCED  AT  GOVERNMENT  EXPENSE, 


Substitute  the  time  derivative  of  unit  vector  to  the  equation  A.20  and  simplify 
then 

R(x.t)  -  (  -  0{  Ry.(x)  +  V  <|>.>'(x)  q^t) }  +  £  <pjx(x)  q^t)  ]  i  (eqn  A.21) 

i  i 

+  [  RX(X)  0  +-  0  ^tPjX(x>  q.(t)  +  V  (p.y(x)  q.(t)  j  j 
i  i 

c.  Kinetic  energy  and  potential  energy 
From  the  equation  A.21 

R(x,t)  =  l  ~  0  (  Ry(x)  +  ^  q.(t)  }  +  V  (p;X(x)  q^t)  ]  i  (eqn  A.22) 

i  i 

+  [  Rx(x)  0  +  0  V  <pix(x>  q;(t)  +  £  (pj>(x)  q;(t)  ]  j 
i  i 

The  dot  product  of  R(x,t)  is 

R(x.t)  •  R(x.t)  -  [  -  0'  Ry(x)  4-  V  <p.>(x)  q.(t) }  (eqn  A.23) 

i 

+  I  <P,X(X)  q.(t)  ]2  +  [  Rx(x)  0  +  9  V<p.*(X)  q.(t)  +  V  v.y(x)  q.(t)  ]2 
1  i  i 

=  02  I  [  Ry2(x)  +-  2  R  (x)  ^  q>iV(x>  qj(t)  +  {  V  (p..v(x)  q.(t)  ] 

»  i  i 

-  2  0  |  {  Ry(x)  +  V  <Pjy(x)  qj(t)  }  X  *jX(x)  q.(t)  ]  +  (I  <PiX(x)  q.(t)  j2 
>  i  i 

+  02  [  Rx2(x)  +  2  Rx(x)  X  (p,x(x)  q.(t)  +-  {  v  q>.x(x)  q((t)  j  2  ] 

i  i 

+  2  0  [  (  RX(X)  +  V  (Pix(x>  q-( t )  }  X  <Pj( X )  q.(t)  j  +  [  V  <p.'(x)  q((t»  |2 

i  i  7  ' 

“  0"  [  !Rx‘(x)  +•  Ry.fc(x>  )  +  2  ^  q,(t)  ;  Rx(\)  <ptx(\i  Rv<\)  (pNxt  ) 

i 

+  1  £  <ptx(x)  q.(t)  }2  +  !  V  <p.>(x)  qjd)  )2  ) 
i  l 
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REPRODUCED  AT  GOVERNMENT  EXPENSE,’ 


+  2  e  V  qj(t)  q;{t)  [  (  Rx(x)  +  <p;X(x)  }  <p;>(x)  -  {  Ry(x)  +  <p;-(x)  )  <p;x(x)  ] 
i 

+  (  (1  <Px(x)  qj(t)  }2  +  {  1  <fyv(x)  q;(t)  }2  ] 
i  i 

From  the  geometry  of  the  system 

Rx2(x)  +  Rv2(x)  =  Rq2(x)  (eqnA.24) 

Rx2(£)  +  Ry:(C)  =  Rfl2(e> 

Substituting  equation  A. 24  into  equation  A. 23  and  apply  the  equation  A. 9 
then  kinetic  energy  is 

I  •  c  t 

T  =  —  02  [  fQ  R02(x)  dm  +  M  R02{£)  +  J0(  (I  <PjX(x)  q,(t)  )2  (eqn  A.25) 

+  (  Z  <Pf'(x)  qs(t)  )2  }  dm  +  M  {  (  Vtp^C)  q.(t)  )2 

i  i 

+  (  I  <Piy(*)  qj(t)  )2  } 

i 

C 

if  +  2  1  q/t)  (  J  (  Rx(x)  9jX(x)  +  Ry(x)  9jy(x)  )  dm 

% 

+  M  (  RX(C)<P;X(C)  +  Ry(C>  <PiV(C) }  1 

1  l 

+  r  [f  ((I  <P;X(x)  qj(t)  )2  +  (  V  9>(x>  q.(t)  )2  }  dm 

2  U  j  • 

+  M  {  (  v  <px(C)  q^t)  )2  +  (  V  9j>‘(C)  q^t)  )2  )  }  ] 

i  i 

•  .  ,  e 

+  6  I!  q^f)  (  J0»  Rx(x»  <pi>(x)  -  Ry( x )  tp^x)  }  dm 

+  M  (  RX(C)  <p,y(C)  -  R>(C)  <PjX< C)  } 

C 

+  1  I  J(|(  <PjX(x)  ~  9j'(x)  <PjX(C>  dm 
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+  M  (  <pjX(C)  <p;y(0  -  <Pj>(C)  )  }  ] 

1  •  •) 

+  - Ir  02 

2  rzz 

Now  let's  apply  the  orthogonality  relationship 

£ 

J  {  <PjX(x)  <PjX(x)  +  <p.v(\)  <p.-v(x) }  dm 

+  M  (  (p^C)  <pjx(£)+  <Piy(  C  >  <pjv(€) }  -  0  (for  i  *  j  ) 

=  Mj  (for  i  =  j  ) 

to  the  equation  A. 26  then 

T  -  ~T  ®2  I  j!  ro2(x>  dm  +  M  fy2**)  +  I  V«)  Mi 

“  u  i 

£ 

+  2  Z  ^i(l)  {  f0<  Rx(x)  <PiX(x)  +  Ry(x)  q>iy(x) )  dm 

+  2M(  RX(C)  <PjX(C)  +  Ry(C)  <pjy(  C) ) }  +  I  ) 

+  4-iv«> 

“  i 

•  •  E 

+  0  Z  I  JQ(  Rx(x)  <Piy(x)  -  Ry(x)  <PiX(x) }  dm 

+  M  {  Rx(C)«Pi-v(t)  -  Ry(C)  <Pjx(C) } 

£ 

+  1  (  J  <  <PjX<x)  <pj-v(x)  -  «Pjy(x)  <PiX(£)  )  Am 

■  0 

+  M  ( <p.x(C)  <p,y(C)  -  <p.>(C)  <p,x(C) ) }  ] 

From  the  equation  2.25,  potential  energy  is 


(eqn  A.  26) 


(eqn  A. 27) 
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(eqn  A.  15) 


u  =  — t~Z  «i2  Mi  qfd) 

-  i 

d.  Lagrange's  equation 

Apply  Lagrange's  equations  A.l  and  A. 2 

then 

3T 

— 0 

-75-  -  ®  I  /  R0:'x> dm  +  m  R0:<f)  +  I  q,2(t)  m, 

_  ,  f 

+  2  Z  q^)  (  J  J  Rx(x)  <PiX(x)  +  Ry(x)  <Piy(x)  )  dm 
i  u 

+  2  M  (  Rx(C)  <pA<£)  +  Rv(£)  <p>(£) ) }  +  1  ] 

ZZ 

£ 

+  I  qj(0  1  JQ  (  Rx(*)  <Piy(x)  -  Ry(x)  (p^x) }  dm 
+  M  {  RX(C)  <prv(C)  -  Ry(  £>  <pix(C) } 

+  Z  q/1)  (  J0(  <PjX(*>  <?/(*)  -  tp^x)  (pj^f)  )  dm 

+  M  (  ip.\l)  <Pj-v(£)  -  <pjy(  C)  <p;x(C) ) }  ] 

—  t  *7$-  1  "  0  I Jjj  R02(M dm  +  m  Ro2^)  +  I  q|2(t)  Mj 

+  2  Z  q/O  {  J0  (  Rx(x)  q)ix(x)  +  Ry(x)  <pjy(x)  )  dm 
+  2  M  (  Rx(C)  <pix(C)  +  R  (C)  q>.-v(C)  )  }  +  I.  j 

*  ZZ 

+  2  e  Z  qj(t)  q;(t)  Mj 
i 

•  .  C 

+  2  0  Z  q;!  J()(  Rx(x»  <PjX(x)  +  Ry(x)  (p(-v(x)  )  dm 
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b 


V 


+  I  <£(0  I  fg  {  RX(X)  %'  (X)  ~  Ry(X)  <PiX(X)  }  dm 

+  M  (  RX(C)  <P'y(C)  -  Ry(C)  <pjx(  £) } 

l 

+  X  qj(t)  ‘  f0<  <pjX(x)  Vf(x)  ~  <Pjy(x)  <PiX(£) )  dm 
+  M  (  (pjX(C)  <Pj-v(C)  -  cpj^C)  <PiX(C) ) }  ] 

dU  cl  d  <?U 
-  =  — _  =  - r  — ]  =  0 

cQ  c0  dt  cO 

— r—  =  02  [  qh(t)  Mh  +  {  J0  (  RX(X)  <pix(x)  +  Ry(x)  q>.>'(x) )  dm 
+  M  {  R  (C)  <p-x(£)  +  R,.(C)  <p.y(£) ) }  +  I,  ] 

3  1  zz 

+  0  1  5i(t)  [  <ph*(x)  <Piy(x)  -  <ph'v(x)  tp.x(x)  dm 

+  M  ( <phx(C)  <prv(£)  -  <ph-v(C)  <p.x(C) }  ] 

dT  •  J 

— rr~  =  0  (  i  Q{  rx(x>  %;(*)  +  Ry(x>  <Phx(x) )  dm 

+  M  (  RX(C)  <ph-v(C)  -  Ry(C)  <ph*(£) } 

£ 

+  I!  q/t>  {  JQ(  (P^x)  <ph>(x)  -  Vj>(x)  (phx(£)  dm 
+  M  { <PjX( f >  <PjV(f)  -  <Pj'(f)  <PjX(f) }  I  +  qh(t)  Mh 

-7-  [  -TS-  1  *  0  [  f  R.(x)  <ph-'(x)  +  R  (XMP,  x(\>  )  dm 

dt  C qh  ‘0  x  n  >  " 

M  (  RX(C>  <phy(C)  -  R  (C>  <phx(C)  } 


s 
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+  I  q/0  '  L(  <PjS(x)  <Phy(x)  -  <pj-(x)  tphx(f)  Jni 

j 

+  M  1  q>.x(C)  <pj>‘(C)  -  <Pjy(C)  <p,x(C)  }  ]  +  qh(t)  Mh 

#  £ 

+  0  J  ^j(t)  l  10  '  %*(*)  “  <Pjy<x)  <phx<x> !  dm 

+  M  [  <p  *(£)  <phy(C)  -  <pjV(£)  <Phx(C)  1  +  <lh(l)  N,h 

au  -> 

—  =  co-h  Mh  qh(t) 
dclh 

4.1.0 

dt 


Plug  these  quantities  in  the  equation  A.l 
then 

0  |  j  R02(x)  dm  +  M  R02(C)  +  Y,  qflt)  M; 
i  {)  i 

f 

+  2  ^  q.(t)  {  Jy(  Rx(x)  (PiX(x)  +  Ry(x)  (Pjy(x)  )  dm 

+  2  M  (  RX(C)  <PiX(C)  +  R  (C)  q>iy(C) ) }  +  1  ] 

*  2  0  y  qL(t)  qj(t)  M; 

i 

? 

+  2  0  ^  qi(t)  !  J  (  Rx(x»  (p^x)  +  Ry(x)  (p^x)  )  dm 

l  } 

+  2M(  R  (£)  9;x(£)  +  fy £)  <P.y(f)  > }  +  Ir  02  1 

A  i  j  *zz 

*  v  q(t)  {  |  *  <p,x(x>  <p.>(x)  -  <Pjy(x)  9,x(£)  }  dm 

+  M  ( (p,x(£»  <Piy(£)  -  <p,v(£>  <Pixi£)  1  =  t 


(eqn  A.2S) 
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and  the  equation  A. 2  becomes 


6  [  J*  (  Rx(x)  <ph>(x)  +  Ry(x)  <phx(x) )  dm  (eqn  A.29) 

+  M(  Rx(t)  <Phy(0  -  Rv(f)  <Phx(f) } 

£ 

+  I  «lj(t)  {  JQ(  <PjX(x)  <Phy(x)  -  9>(x)  <phx(t)  dm 

+  M  {  <PjX(£MPiV(0  -  <Pjy(C)  <PiX(C)  }  J  +  qh(t)  Mh 
• .  f 

“  6“  [  qh(t)  Mh  +  (  J  (  Rx(x)  <phx(x)  +  Ry(x)  <phv(x)  )  dm 
+  M  1  RX(G)  <phx(C>  +  Ry(C>  <ph>{£) ) }  ] 

+  2  0  V  q^t)  [  j^(  <pjX(x)  <ph-v(x)  -  (p>(x)  <phx(x)  }  dm 
+  M  (  <p.s(t)  <phy(£)  ~  <Pjy(£)  H  +  «2h  Mh  qh(t)  =  “  2  C  Wh  Mh  qh(t) 
(  h  =  1,  2,  3 . n  ) 


e.  Results 

Table  5  and  6  show  the  eigenvalues  of  Aluminum  Alloy  and  Composite 
material  boom  from  NASTRAN  simulation  results  and  Figure  A. 3  shows  first  two 
mode  shape  of  the  double  link  boom  in  planar  motion.  As  shown  in  Figure  A.4, 
applied  torque  to  maintain  15  r.p.m  of  rotating  speed  increased  to  one  hundred  times 
of  the  torque  that  applied  to  the  same  boom  in  3-dimensional  motion  because  its  mass 
center  is  far  away  from  rotating  center.  Angular  displacement  change  with  time  was 
almost  same  as  3-dimensional  case.  Displacement  in  x  and  y-direction  was  much  bigger 
than  that  of  3-dimensional  motion  because  tip  position  is  far  away  from  the  rotating 
center  and  concentrated  mass  attached  to  the  tip. 

While  rotating  speed  increases,  boom  deflects  negative  y-direction  and  positive 
x-direction  until  reaching  its  equilibrium  position.  After  cutting  ofT  torque,  boom 
remains  equilibrium  position  and  oscillates  with  some  amplitude.  Magnitude  of 
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NSE, 


deflection  at  equilibrium  position  was  about  6.2  inchs  and  slope  at  the  tip  position  was 
2.3  °.  These  values  are  much  bigger  than  that  of  3-dimensional  motion  but  it  is  quite 
natural  because  of  its  big  radius  of  rotation. 

Figures  A. 5  -  A.  12  show  angular  displacement,  angular  velocity,  generalized 
displacement,  deflection  in  each  direction,  magnitude  of  deflection  and  slope  at  the  tip 
position  respectively. 

TABLE  5 

REAL  EIGENVALUES  OF  ALUMINUM  ALLOY  (2D) 


mode 

no. 


radians 

o>i 

cvcles 

o>i 

eeneralized 
w  mass(  \L) 

generalized 

stillhess(Kj) 

3.56SS16E  +  00 

5.679947E-0I 

1.000000E  +  00 

1.273645E  +  01 

1.796733E  +  01 

2.8595S9E  +  00 

_ 

1 . 000000 E  +  00 

3.22S24SE  +  02 

TABLE  6 

REAL  EIGENVALUES  OF  COMPOSITE  MATERIAL  (2D) 

radians 

"i 

cvcles 

«i 

generalized 

mass(Mj) 

eeneralized 
stifl'nessi  Kj) 

4.46132SE  +  00 

7.100423E-01 

1. 000000 E  +  00 

1.990344E  +  01 

2.2461 15E  + 01 

3.574S04E  +  00 

1 .000000 E  +  00 

5.045035E  +  02 
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Figure  A. 4  Applied  torque  vs.  time. 
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APPENDIX  B 

DETA.LEDDER.VATIgM^WlvSgfl^UATIONS 


FOR  THREE 


Apply  Lagrange's  equations  2.7  and  2.S  into  equations  2.28  and  2.25 

then 


dT 

cQ 


=  0 


=  0  [  V  q.2(t)  M.  +  f  V2(x)  dm  +  M  Rx(£) 
i  u 

e 

+  2  I  q/t)  (  J  Rx(x)  <P;x(x)  dm  +  M  Rx(£)  <PiX(C>  } 


-  2  (  J  \  £  <p;z(x)  q.(t)  )2  +  M  (  Z  (Piz(£)  q.(t)  )2  +  Ir  ] 
U  i  •  zz 


+  £  I  ^(0  q/0  (  JlV,x(x)  9jy(x)  dm  +  M  <pjx(£)  9jy(£) } 

i  j 

+  I  qs(t)  {  J0\(x)  <Piy(x)  dm  +  M  Rx(£)  <prv(£) } 


]  -  e  [  y  q2(t)  M.  4-  jVx2(x)  dm  +  M  Rx2(£) 

£ 

+  2  £  q^t)  {  JQ  Rx(x)  <pjX(x)  dm  +  M  Rx(£)  <psx(  £ >  ) 

~  2  {  J  J  I  <PjZ(x)  qj(t)  )2  +  M  {  I  %z(l)  q;(t)  )2  +  1^  ] 

+  2  0  £  qj(t)  [  qj(t)  M.  +  Jq£Rx(x)  y*(x)  dm  +  M  Rx(£>  tp^f)  } 

~  2X q/1)  ( J£<Piz(x) <Pjz(x> dni - M  <Piz(f) q>j7(f )  >  ] 

••  c 

-  y  q;(t)  [  y  qj(t)  {  Jo  <PiX(x)  <prv(x)  dm  +  M  <PjX(£)  <PjX(f) 


120 


\ 

) 


£ 

-  J  <piy(\)  <Pj'x(\)  dm  -  M  <Pj-v(C)  <p-x(  C ) 

£ 

-  _[  Rx(x)  cpj- (x>  dm  -  M  R^fHp^D  ] 

f 

=  e-  [  qjtn  Mj  +  J  RX(X)  <PjX(x )  dm  +  MRX(C)  <p*(£) 

£ 

“  2  1  J  <PiZ(x>  (phx(X!  dm  +  M  <p;z(f)  <p,/(f)  }] 
i  '* 

»  »  £ 

*  6  1  Mjlt)  [  J  Jphx(\)  <pjV(x)  dm  +  M  <phx(fi  <p.y(£) 

£ 

-  J^(pi>(\)  <p-x(x)  dm  -  M  <pj-v(C>  (p.x(E^  ] 


=  qh(t>  Mh 

*  •  £ 

*61  q,(t)  [  J^<pix( x >  (ph>(x>  dm  +  M  y  *(t)  <ph-v(C) 

£ 

-  f  «p.-'(x>  (phx(x)  dm  -  M  <p>(£)  <phx(£)  ] 

•  £ 

*  0  [  J[  Rx(x)  <phy(x)  dm  +  \1  Rx(f)  <phy(£)  ] 


=  Mh(t)Mh 


..  £ 

+  9  [  1  qj(t)  {  J(<pjXixi  (phy(x)  dm  +  M  <pix( £ >  <ph>(£) 

£ 

-  J  <Pj'(x)  <phx(x>  dm  -  M  <p.y<f>  <phx(f)  } 
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£ 

+  J0Rx(x)  (phy(\)  dm  +  M  Rx(f)  9h>'(«)  ] 

•  .  C 

+  0  Z  ^(0  I  Jn<piX(x)  Vhy(x)  dm  +  M  <Pj*(C)  <ph'(C) 
i  u 

£ 

-  J0<Piy(x)  <Phx(x>  dm  -  M  <prv(f)  (phx(f)  ] 


cU  _  dU  _  d  f  cL' 

To  TS  dt-1  1  =  ° 


au 

e% 


<02h  Mh  qh(t) 


d\J 

£% 


d  dV 
— (tt-i  =  0 

dt  ^qh 


Now  plug  all  these  quantities  in  equations  2.7  and  2.8 
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APPENDIX  C 

NASTRAN  PROGRAM  FOR  DYNAMIC  (  MODAL  )  ANALYSIS 


1.  DYNAMIC  ANALYSIS  IN  PLANAR  MOTION 


S  Dojble  link  flexible  boom  for  dynamic  analysis!  14  grid  points ) 

S  in  2  dimension  space  with  tip  mass  (37.5  lb.  )  linkl :  14  ft,link2:12  ft 
S  Material: aluminum  allov 
S 

^  *  *  *  *  *  *  *  *  *  *******************  *  *  *  * 

s 

S  EXECUTIVE  CONTROL  DECK 
S 

c  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  4c  *  *  *  *  *  *  *  *  *  sjt  *  *  *  *  * 


id  kang.dvnamics  SN-ROSS 

sol  3  ~  ' 

time  10 

diae  8.13 

cent! 

S 

^  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  4c  * 

S 

S  CASE  CONTROL  DECK 

S 

^  **  ****  .■*$**  .*  **. It  *«***. ■**!!!. *******  4 *  * 

S 

title  =  Modal  analysis  of  double  link  flexible 

title  =  boom  in  planar  motionf  Aluminum  allov) 

echo =  both 

method  =  120 

snc  =  101 

disp  =  all 

output!  plot) 

plotter  nastran 

axes  z.x.v 

cscale=  1.8 

view  O..0..O. 

paper  size  14.0  bv  10.0 

set  1  =  all 

find  scale,  origin  l.set  1 

plot  modal  delormation  0,set  1, origin  1, shape 

maximum  deformation  5 

begin  bulk 

S  “ 

^  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  4c  *  *  *  *  *  *  *  * 

S  BULK  DATA  DECK 

S 

^  &  *  *  **********  *  *  *  $  4c  ********  *  *  *  * Jl«  ♦ 

S 

S  Define  new  coordinate  system  for  convenience 

Scord2c.l..0..0..0.,0.,0.,l..  +  23 
S  +  23.1..U...O 

S . 

S  grid(node)  data 

S--- . 

grid, l.jO. .0. .(N  t 


S  element  data 

S . 

cbar,2l.l02.1.2.0..0.,l. 

).-,*( iV( l - 

S  Change  these  boom  properties  when  material 
S  changes  AL.  to  COMPOSITE 
8 

S  —  DATA  FOR  ALUMINUM  ALLOY  — 

pbar.  102, 103.1. 2101.1. 2447, 1.2447, 2.4894 

matl.l0ll.0\  +  7.3.}  +  6„2.5362-4 

conm2.103.14..9.7176-2 

spcl.1 1 1.1^3456,1 

spcl.l I2.345.2,tfrru,14 

spcadd.101.1 1  f.l  12 

eigr,120,meiv„„7 

Sparam.autospc.yes 

S  —  DATA  FOR  COMPOSITE 
S 

pbar.  102, 103.2.01 68. 1 .945  U  .945 1 ,3.8902 

mat  1.103, 1.01  +  7  .25.2.5362-4 

conm2.103.14..9.'/\76-2 

spcl.l  11.123456.1 

spcl.l  12.345. 2,thru.l4 

spcadd.101,1 1  f,l  12 

eigr.l20.meiv.,„7 

Sparam.autospc.ves 

enddata 


2.  DYNAMIC  ANALYSIS  IN  3  DIMENSION  SPACE 


S  Double  link  flexible  boom  for  dynamic  analvsis(  14  grid  points ) 
S  in  3  dimension  space  with  tip  mass  (37.5  lb'  ) 

S  link  1:14  ft.link2:12  ft 

S  anale  between  LI  and  x-axes  is  70, LI  and  L2  is  126  degree 
S  Material:  Aluminum  allov 
S 

^444$* 4 4 4 4 4 4 4  4 4  4  4 4  4  4  4  4  4  4  4  4  4  4  4  4  4  4  4  4  4  4 

S 

S  EXECUTIVE  CONTROL  DECK 

s 

I *  *  *  *  * *  *  *  *  *  *  *  * * ******  *  *  *  *  *  *  *  *  *  *  *  *  *  *  * 

S 

id  kang, dynamics  Snross 

sol  3 

time  10 

diaa  8,13 

cerid 

c *************** ******************** 

S 

S  CASE  CONTROL  DECK 
S 

^  ******  * *  *  *  *  *  *  *  *  * *  *  *  *  *  *  *  *  *  ********** 

s 

title  =  2  link  flexible  boom  in  3  dimension! A1.COM) 

echo  =  both 

method  =  120 

spe  =  101 

disp=  all 

output(plot) 

plotter  nastran 

cscale=  1.8 

view  0..0..0. 

paper  size  14.0  bv  10.0 

set  1  =  all 

find  scale,  oriain  l.set  1 
maximum  deformation  10 
axes  my.x.z 

plot  modal  deformation  O.set  1. oriain  I. shape 
axes  x.y.z 

plot  modal  deformation  O.set  1. oriain  l, shape 
becin  bulk 

^  4.  4  4  4  4  -4  -4  4-  4  4  4  4  4  4  4  4  4  4  4  4  4  4  4  4  4  4  4  4  4  4  4  4 

S  BULK  DATA  DECK 
S 

^  4  4  4  4  4  4  4  4  4  4  4  4  4  4  4  4  4  4  4  4  4  4  4  4  4  4  4  4  4  4  4  4 

S 

S  Define  new  coordinate  system  for  convinience 
S 


Scord2c.l..0..0.0.,0.,0.,l.,+  23 
S  +  23,.5,0..  S..67 


S  grid( node)  data 
S 

arid. I,  .0..0..0. 

=  .*(  1 ).  =  .-(S.2085),  =  ,*22.5526 

arid. 9.  .44.O3SS.0..  177.7653 
=  /:( 1  >.  =  .*(- 1 3.4*06).  =  ,*  19.8969 

S 

S . 

S  element  data 


cbar,2 1.102.1 .2.0..0. ,  1 . 

^  ^1  ),“.*(  1  ).*(  *).““ 

s 

S  Change  following  datas  when  the  material 
S  changes  from  AC.  to  COM. 

S 

S  —  DATA  FOR  ALUMINUM  ALLOY  --- 
S 

pbar.  102. 103. 1 .2 101 . 1 .2447, 1 .2447.2.4894 
matl. 103.1.01  -+-  7  3.)  +  6, .2. 5362-4 
conm2.103.14.j9./ 176-2 
sjsc  1.1 01, 1234o6.1 
Ssncl.l  1^.345.2,thru.l3 
Sspcadd.101.11  (.112 

eigr.120.mgiv..., 7 
Sparam.autospc.yes 

S  —  DATA  FOR  COMPOSITE  — 

nbar.  102. 103. 1 .2 101 .1 .2447. 1 .2447.2.4894 

matl. 103. 1.01  +  7.3. t  +  6, .2. 5362-4 

conm2. 103. 14tJ9.7 1 76-2 

SDcl.loi.l234o6.i 

Sspc  1.1 12. 345.2. thru,  13 

Sspcadd.101.11 1.1 12 

eier.120.mgiv..  ..7 

Sparam.autospc.yes 

enddata 


APPENDIX  D 

DSL  PROGRAM  SOLVING  THE  DYNAMIC  EQUATIONS  OF  MOTION 
1.  PROGRAM  OF  DOUBLE  LINK  FLEXIBLE  BOOM  IN  PLANAR  MOTION 


*  This  program  solves  the  dynamics  of  a  dobule  link  flexible  boom  system 

*  in  planar  motion.  The  applied  torque  (  t )  was  changed  to  see  the  effects 

*  of  the  torque  and  also  damping  coefficient  (  £  )  was  varied  0.0  to  0.5  %. 

*  This  program  automatically  calculates  the  slope  and  deflection  of  the  tip 

*  position  in  each  direction  (  local  x,  y,  z  )  from  simulation  results. 

*  _ _ _ _ _ _ _ _ _ _ 

*  THE  FOLLOWING  PARAMETERS  ARE  DEFINED 


N  ;  number  of  discritized  boom  element 

FL1  ;  length  of  lower  boom 

FL2  ;  length  of  uppper  boom 

TM  ;  tip  mass 

RHO  ;  mass  per  unit  length 

DX  ;  differential  length  between  each  grid  point 

RX  ;  local  x  distance  from  origin 

TAO  ;  applied  torque 

OMG1  ;  1st  mode  natural  frequency 

OMG2  ;  2nd  mode  natural  frequency 

PX1  ;  1st  mode  shape-displacement  in  local  x  direction 

PY1  ;  1st  mode  shape-displacement  in  local  y  direction 

PZ1  ;  1st  mode  shape-displacement  in  local  z  direction 

PX2  ;  2nd  mode  shape-displacement  in  local  x  direction 

PY2  ;  2nd  mode  shape-displacement  in  local  y  direction 

PZ2  ;  2nd  mode  shape-displacement  in  local  z  direction 

ALP  ;  angle  between  upper  boom  and  local  x-axes 

RX1L  ;  1st  mode  shape-rotation  at  the  tip  position  w.r.t  x-axes 
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RZ1L  ;  1st  mode  shape-rotation  at  the  tip  position  w.r.t  z-axes 

RY2L  ;  2nd  mode  shape-rotation  at  the  tip  position  w.r.t  y-axes 

RZ2L  ;  2nd  mode  shape-rotation  at  the  tip  position 

TH  ;  angular  displacement 

THD  ;  angular  velocity 

Q 1  ;  1st  mode  generalized  coordinatest 

Q1D  ;  time  derivative  of  1st  mode  generalized  coordinates 

Q2  ;  2nd  mode  generalized  coordinates 

Q2D  ;  time  derivative  of  2nd  mode  generalized  coordinates 

ZETA  ;  damping  coefficient 

R1RZZ  ;mass  moment  of  inertia  for  R.F  electronic  box 
R1BZZ  ;mass  moment  of  inertia  for  boom  and  tip  n  mass 
Cl  ;  starting  up  time  coefficient 
C2  ;  magnitude  of  applied  torque 
C3  ;  required  time  to  reach  required  R.P.M 


^  ^  ^  ^  +  -it  +  *  * .+  *  *  *  .^  *  ^  *  «  it  *  *  *  «  +  +  +  *  +  *  *  *  +  *  *  *  +  *  *  ♦  ♦  +  ♦  «  «  +  *  it  *  .-{x  +  *  ii  it  it  ♦  «  *  *  *  it  it 

* 

*  SIMULATIOMS  OF  DOUBLE  LINK  FLEXIBLE  BOOM  IN  PLANAR  MOTION 

* 

*  %  *  *  *  ?;?  *  *  }‘,i  #  #  4  ft  ft  ft  ft  ft  *  ft  ft  ft  ft  *  *  ft  *  ft  *  *  ft  ft  ft  ft  ft  ft  *  *  *  ft  *  *  ft  *  ft  *  ft  ft  *  4c  *  ft  *  *  *  *  *  ft  ft  *  *  *  *  ft  *  *  ft  ft  £  ft  &  «  &  $  ft  *  ft  ft 

4 

TITLE  SOLUTION  OF  SIMULTANEOUS  DIFFERENTIAL  EQUATIONS 
TITLE  FOR  DOUBLE  LINK  FLEXIBLE  BOOM  IN  PLANAR  MOTION 


FIXED  IER.  I  PVT.  N.  I 

CONST  FL1  =  16S.,FL2  =  144..RHO=  3.0690E-04.N  =  13,DELT=  .04 
PA  RAM  Cl  =  40.  ,C2  =  1000.  ,C3=  19.65 


INITIAL 


TM  =  9.7176E-02 
ALP  =  54. 

RIRZZ=  31.0559 
RZ1L=  1.527773E-02 
RZ2L  =  -3. 149482E-02 


BY  CHANGING  THE  ZETA  VALUE  THE  DAMPING  FORCE 
WILL  BE  CHANGED 


ZETA  =  0. 

ZETA  =  0.002 
ZETA  =  0.005 


DIMENSION  SIZE  SHOULD  BE  EXPRESSED  BY  NUMBER 
INSTEAD  OF  CHARACTER 


D  DIMENSION  RX(  13).RY(  13).P.X1(  I3),PX2(  13).PY1(  13),PY2(  13 >.A(3,31 
ARRAY  IPVT(3).B(3) 
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■'  «/■ 


*  REMOVE  THE  ASTRICKS  IN  DATA  STATEMENT  FOR 

*  EXECUTION  OF  EACH  MATERIAL 

j{; 

D  DATA""rX  *24  *48:;7Y*96:*i20:;i^7l68"l  8"2."i668*,  196.2  MfilOjiT;" 

D  ft224. 428. 238. 535, 232. 642 

D  DATA  RY  7*0.0, 19.4164, 3S.S32S, 58.2492, 77.6656, 97.082,1 16.498: 
DATA  FOR  ALUMINUM  ALLOY  =  =  =  =  =  = 


D  DATA  PX1  -4.922269E-06, 

D  ft  -9.844536E-06. 

D  a  -1.4766SOE-05. 

D  =  -1.96S907E-05. 

D  ft  -2.461 132E-05. 

D  =  -2.95335SE-05, 

D  =  -3.4455S3E-05, 

D  ft  -2.457687E-01. 

D  ft  -5.0S7S06E-01, 

D  ft  -7.853463E-01. 

D  ft  - 1. 071851E  + 00,-1. 364817E  + 00,-1. 660928E  + 00' 

D  DATA  PY1  2.800838^-02, 

D  ft  1.089855E-01, 

D  ft  2.383627E-01, 

D  ft  4.1 15835E-01, 

D  ft  6.241 170E-01, 

D  tt  8.714750E-01, 

D  a  1.149234E  +  00, 

D  ft  1.327774E  +  U0, 

D  ft  1.518872E  +  O0, 

D  ft  1.719S15E  +  00, 

D  »  1.92797SE  +  00,2. 140S36E  +  00.2.355983E  +  00/ 

D  DATA  PX2  1.657409E-04, 

D  «  3.314SIOE-04. 

D  «  4.972196E-04, 

D  =  6.629559E-04. 

D  a  S.286S90E-04. 

I)  a  9.944 183E-04, 

D  «  1.160143E-03, 

D  =  4.385572E-02. 

D  ft  2.710247E-01. 

D  a  6.4S7206E-01. 

D  ft  1. 14045 7 E  +  00. 1.708 1 S0E  +  00,2.3 13523E  +  00.' 

D  DATA  PY2  1 .42543 1 E-01 . 

D  ft  5.101371  E-0 1 . 

D  =  1.013169E  +  00, 

D  =  1.563418E  +  00, 

D  =  2.075334E  +  00, 

I)  =  2.467438E  +  00, 

D  =  2.663624E  +  00. 

D  ft  2.632S37E  +  O0. 

D  ft  2.46S005E  +  00. 

D  =  2.103798E  +  00, 

D  ft  1.836722E +  00,1. 424425b +  00.9.847735E-01 

D  DATA  OMG1  3.5688 16E  +  00  ,OMG2  1.796733E  +  01 


=  DATA  FOR  COMPOSITE  MATERIAL  =  =  =  =  = 


DATA  PX 1 


-4.615325E-06. 
■9.23065  IF -06. 

■  1.3  8459  7L-05, 
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* 

* 

* 

* 

* 

* 


u 

3 


n 

d 


-1.8461 29E-05, 

-2.30766 1 E-05, 

-2.769 193E-05, 

-3.230723E-05, 

-2.457675E-01, 

-5.087805E-01, 

-7.853472E-01, 

- 1 .07 1 853E  +  00,- 1 . 3648  20 E  +  00,- 1 .660932E  +  00/ 


* 

DATA  PY1, 

2.800845E-02 

* 

a 

1.089857 

E-01 

* 

a 

2.383633 

E-01 

* 

T* 

4.115847 

E-01 

* 

jj 

6.241187 

E-01 

■/. 

... 

d 

S. 714775 

E-01 

* 

1.149238 

E  +  00 

4 

1.32777S 

E  +  00 

2£ 

1.518877 

E  +  00 

* 

d 

1.719S20 

E  +  00 

1 .927984E  +  00,2. 140842E  +  00.2.355990E  +  00/ 


DATA  PX2,  1.554086E-04. 
a  3.10S165E-04, 

=  4.662230E-04, 

a  6.216275E-04, 

a  7.770293E-04, 

a  9.324276E-04 

a  1.087822E-03, 

a  4.377849E-02, 

a  2.709504E-01, 

a  6.48656  IE-01, 

a  1 . 140408E  +  00, 1 .708 149E  +  00,2.3 1 35 14E  +  00/ 


DATA  PY2/ 


j. 

*  = 

*  2* 

*  ^ 


1.425488E-01 
5.  IO1575E-01 
1.OI3209E  +  00 
1.563480E  +  00 
2.07541 6E  4- 00 
2.467535E  +  00 
2.66372SE  +  00 
2.632930E  +  00 
2.46808  2  E  +  00 
2. 193S55E  +  00 

1.836756E  +  00.1.424434E  +  00.9.847565E-01 


D  DATA  OMG1  4.46132SE  +  00.  .O.MG2  2.2461 15E  +  01 


THESE  CALL  STATEMENTS  CALCULATE  ALL  CONSTANT  TERMS 
AND  COEFFICIENT  INVOLVED  IN  THIS  PROGRAM. 


CALL  FIBZZ(RHO. 

CALL  ONE(RHO,R.  _  .  _  .  _ 

CALL  T H R I:  E(N?fOC  I R\\ f*>U JP^Yl ?PX  2 .  PY2 .TM.S13.Sr4.ST7.ST8. Si  ll. 

Stl2.^Tl(>.STISl 
CALL  FOLR(N.RHO.DX.RX.RY.PXI.PYl.PX2.PY2.SK\Sr9> 

CALL  FIVECN.RHO.DXJRX.RY.PXl.PYl.PX2.PY2.ST6^1  10  ,  „ 

CALL  SIX(N,RHO,DX,fM.RX.ftY.PXl,PVl.P\2.PY2.Sr 19.S  1 20.S  1 21.S  T22) 

DERIVATIVE 

NOSORT 

Tl  =  5. 

TAOl  =  Cl  *TIME**2 
TA02-C2*M--  STEP(C3  )) 

TAO  =  SW ITCIK  I  I M  E.  L  E.T I  .TAOl  .TA02 )  , 

A(  1 .1  )=  R1 BZZ  +  RI RZZ  +  X2S'*2  +  X3**2  +  (ST1  +  S I  3)*X2  +  (S  T2  +  S 1 4)*\3 
A  1.2)  =  SI  5  +  ST7  +  ( S  16  +  S  TS)*X3 
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A(  1.3 
At  2.1 
AP  ? 


A 

A 

A 

A 

B< 

B 


2.j 


3.1 

■> 


j.j 


=  ST9  +  ST1 1  +JTST10+  ST12)*X2 
■  =  ST5  +  ST7  +  (ST6  +  STS  )*X3 
=  1.0 
I  =  0.0 

i=  ST9  +  ST1 1  +(ST10+  ST12)*X2 

i  =  0.0 
1=  1.0 


1 )  =  TA0-X1D*(X2D*(2.*X2  +  ST1  +  ST3)  +  X3D*Q.*X3  +  ST2  +  ST4)) 
2  =  X 1  D**2*(X2  +  ST20  +  ST16)-2.*X  1D~(X3D*(ST19  +  STS))-OMG  1* 
-2.  *ZETA*OM  G 1  *  X  2  D.  ’  _ 


B(3)  =  X1D**2*(X3  +  ST22  +  STIS)-2.*XID*(X2D*(ST21 
-2.*ZETA*OMG2*X3D 
CALL  DGEFA  (A.  3.  3.  IPVT,  IER) 

IFUER.XE.Oj  GOTO  112 

-  -  -  lpvT  B 


STl2))-OMG2**2*X 


0) 


*120 


CALL  DGESL  (A.  3.  3. 

X1D=  1XTGRL(0.,B(1 
X2D=IXTGRL  0..B  2 
X3D=1NTGRL  0..B13 
XI  =  1NTGRL(0..X1D) 

X2=IXTGRL  0..X2D) 

X3-INTGRL  0..X3D) 

TH  =  XI 
THD=  X1D 
"  1  =  X2 
1D  =  X2D 
2=  X3 
2D  =  X3D 

VVX=PX1(X)*Q1  +  PX2(X)*Q2 
WV=  PYlfX'WSl  +  PY2fX)*Q2 
WXY=  \VX  **2  +  WY**2 
W=  SQRTLWXY) 

SLOP  =  (RZlL*Ql  +  RZ2L*Q2)  57.2957 
WRITE(6.120)RIBZZ.STl.ST2.Sl  3.ST4.S  17  ST8.ST1 1. ST  12. ST  16. ST  18. 
ST5.ST9.ST10.STi9.Sf20.ST21.ST2i 
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FORMAT!  2X,1S(F12.5.2X)) 
RFTL’RN 

12  WRITE(6.i  14)  TIME. IER 
14  FORMaY( 2X.  0  IER  =  .17) 
CALL  EXDJOB 

WR  Ij'F.l  6.J  1 6 )(QU  F 1 ).  I F 1  =  1  .X ) 
2\..' FS.3.2X) 


*116  FORMAT! ' 


PRIXT  TILTH  D.  Q1.Q2AYX.WYAV.SLOP 
COXTRL  YlXTI\F=  40..DELPRf  =  .20 
SAVE  0.025.TH.THD.Q1.Q2,WX.\VY.W,SLOP.TAO 
HXD 

PA  RAM  C 1  =  SO.  ,C2  =  2000.  ,C3  =  1 1 .50 
EXD 

PA  RAM  Cl  =  160.  ,C2  =  4000.  .C3  =  7.4 

GRAPH  (G10  XT  =  7,LO=  -S00..DE  =  TEK61S.SC  =  800.)  TI.ME(XT  -  5.LE  =  10.  .L  X  =  ... 
SEC  ),TAO(LX  =  LBS-IX  .RL  =  l.LI=  1)  ... 

.TAOf  1*0  =  0..AX  =  OMIT.RU  =  2. LI  =  2)  ... 

.TAO(PO  =  0..AX  =  OMIT.RL  =  3. LI  =  3) 

GRAPH  (Gl.XI  ®  7.LO  =  0..DE  =  f EK61S,SC  =  8.  )  TIME(XI  ®  5.LE*  10..UX-  SEC  )  . 


,TH(LX  =  RAD.RL  ®  l.LI  =  1)... 

,TH(PO®0.,AX  =  OMIT.RL  =  2. LI  =  2) ... 

,TH(PO  =  0.,AX  =  OMIT.RL  =  3. LI  -  3) 

GRAPH  (G2.XT  =  7.LO  =  0..DE  =  TEK61S.SC  =  .3  )  TIME(XT  =  5.LE  =  10..L X  =  SEC  ) 
,THt>(LX=  RAD  SEC. RL=  l.LI  =  1)  ... 

,THD<  PO  =  0..AX  =  OMiT.RL  =  2.LI  =  2)  ... 

.THD!PO  =  0..AX  =  OMIT.RL®  3.LI  =  3) 

GRAPH  <G3.XT  =  8.LO®  -8.0  .SC®  1.0  .  DE  =  TEK6I8)  ... 

TIME!  XT  =  5.LE=  10.XX  =  SEC  ).Q1(LX=  TX  .RL  =  1,LI=1),  ... 

Ql(PO  =  0..AX  =  OMlf .RL’  =  2. LI  =  2),  ... 


Ql(PO  =  0..AX  =  OMIT.RL  =  3. LI  =  3) 

GRAPH  (G4. XT  =  SXO®  0.0  .SC  =  .075 

TIMF.t XI  ®  5.1.E  =  10. XX  =  SEO.Q2(LX  =  TX  . 
)2(PO  =  0.,AX  =  OMlf  .RL  =  2. LI  =  2),  ... 

>2(PO  =  <>.  A  X  =  OM I T.  RL  =  3  X I  =  3) 

}4,XT  =  8X0  =  -I.2E-3.SC  =  3.0E-04, 


,DE=TEK618)  ... 

'■  •*'"  RL®  1X1=  i: 


*GRAPIl  (< 


DE  =  TEK618) 
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*  TIME(NT  =  9. UN*  SEC  .  ^ 

GRAPH  (G5.DE  =  TEK61S,L0  =  -0.O.SC*  1.5  ,X 
,LE=  10.).WX(UN='IX.RU  =  l.LI  =  1 ) .. 
WX(PO  =  O..AX  =  OMlf.RU  =  2.LI  =  2) . 
,WX(PO  =  O..AX  =  OMIT.RU  =  3. LI  =  3) 
GRAPH  (G6.DE  =  TEK61S.LO  =  -16.. SC  =  2.0  .NT 
,LE=  10.), WY(  UN  =  TX'.RU  =  l.LI  =  1 )  .. 
WY(PO  =  O..AX  =  OMIT.RU  =  2. LI  =  2)  . 
AVYfPO  =  O.AX  =  OM I  T,RL  =  3, LI  =  3 
*GRAPti  (G7.DE  =  tEK618.LO  =  -.5,SC  =  .13  ) 

y.  i  a'  l  'iv  rrr"  \'i  o  \ 


;\=  in  ) 

0_.SC  =  1.5  ,NT  =  8  )  TIME(NI  =  5.LN=  SEC  ... 


=  8  )  TIME(NT  =  5,UX=  SEC' ... 


)  TIME(NT  =  5,UX  =  SEC'.LE  =  10.)  ... 

*  L’N  =  IN  SEC*  NI  =  8) 

GRAPH  (GS.DE  =  TI~  K618.LO  =  0.0.SC  =  2.5.NI  =  8  )TIME(NI  =  5,UN  =  SEC'  ... 

.LE  =  10.).W(  UN  =  IV. RU  =  l.LI™  0  ... 

AVl PO  =  ()..AX  =  OMIT.RU  =  lLI  =  2)  ... 

.W(PO  =  0.JAX  =  OMIT.RU  =  3.LI=  3) 

GRAPH  (G9.DE  =  1'EK618.LO=  -6.4.SC=  .8. NT  =  8  )  TIME(NI  =  5. UN  =  SEC'  ... 

.LE=  10.).SLOP(UN  =  DEG  .RU  =  l.LI  =  1) ... 

.SLOPlPO  =  0..AX  =  OM IT.RU  =  2. LI  =  2) ... 

.SLOP(PO  =  0..AX  =  OM  IT.RU  =  3, LI  =  3 
LABEL  (G10)  APPLIED  TORQUE 
LABEL  (Gl)  ANGULAR  DISPLACEMENT 
LABEL  G2  ANGULAR  VELOCITY 
LABEL  G3)  GENERALIZED  DISPLACEMENT  Ol 
LABEL  G4  GENERALIZED  DISPLACEMENT  Q2 
LABEL  G5  DISPLACEMENT  IN  EXTENSION 
LABEL  G6[ DISPLACEMENT  IN  TRANSLATION 
■LABEL  (G7)  TIME  DERIVATIVE  OF  03 

LABEL  (GS)  MAGNITUDE  OF  DEFLECTION  AT  TIP  POSITION 
LABEL  G9  SLOPE  AT  THE  TIP  POSITION 
END 
STOP 

FORTRAN 
C 

Q  X  *  *  *  *  *  *  *  *  *  *  *  *  X  *  4c  XX  XX  *  *  *  *  *  *  **********  *  *  *  **********  *  *  *  *  ********  *  *  * 

SUBROUTINE  FIBZZ(RHO.FLl.FL2JALP.TM.RIBZZ) 

£  *  *  *  *  *  *  *  *  *  *  *  *  *  .•;<  *  *  *  *  *  *  *  *  *  *  *  * *  *  *  *  *  *  *  *  *  *  *  *****  *  *  *  *  *  *  *  *  ******  *  *  *  *  *  * 

c 

IMPLICIT  RE.AL*8(.A-H,0-Z) 

PI  =  ARCOS(-l.O) 

ANG  =  ALP-PI  180. 

D 1  SO  =  ( FL 1  +  (FL2  2.)*COSf ANG))**2  +  (( FL2  2.)*SIN(ANG))**2 
D2SQ  =  ( FL 1  +  FL2*COS(  ANG))**2  +  ( FL2*SlN(ANG)r*2 
BM 1  =  RHO*FLl 
BM2=  RHO-FL2 

RIBZZ  =  BM  1*1  LI **2  3.  +  BM2*FL2**2. 12.  +  BM2*D1SQ+  TM-D2SQ 
WRITEI6.il  ItRIBZZ 
111  FORMAT(2X.RlBZZ=.F12.5) 

P  R I N  T  -  '  *  *  *  *  *  *  *  * *  *  *  *  *  *  *  *  ********  ***********  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  • 

RETURN- 
END 

c 

Q *  *  *  *  *  *  *  *  * * **********  *  ********** *  *  ****** *  *  *  *  *  *  *  *  *  *  *  * *  *  ******** *  * * 

SUBROUTINE  ONE(RHO,RX.RY,PXl.PYl.DX.N.STl ) 

£  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  ffi  X  X  *  *  *  *  *  *  x  x  X  X  X  *  *  X  *  *  *  x  X  X  X  x  X  X  *  *  X  X  X  *  X  X  X  X  X  X 

c 

IMPLICIT  REAL*8{A-H.O-Z) 

DIMENSION  RX(N),RY(N),PX!(N).PYI(N) 

SON  El  =0.0 
SONE2  =  0.0 
DO  10  1=  l.N-1 

OXEt  =  RX(I)-PXKI) 

OXE2=RYnV:!PYim 
SONE 1  =  SONE I  +  ONE  1 
SONE2=  SONE2  +  ONE2 
10  CONTINUE 

IT  I  =  ( SONE  1  4-  RXlN  UPX  1(N)  2.)*DX*RHO*2.0 
T 1 2  =  ( S O N E 2  -t-  R Y ( N  r P Y 1  ( X i  2.)*DX*RHO*2.0 

ST1-TI1+  112 
WRITE!  6.1 1 2  >ST  1 


)  TIME(NI  =  5,UN=  SEC'  ... 


132 


no  nn  no  on  on 


112  FORM AT(2X.'ST1  =  '.FI 2.5} 

PRINT  *  *  * *  * * *  * * * * * * ***&************************** * *************  *  *  * 

RETURN 

END 

*  *  *  *  *  *  *  *  *  ft  *  -H  *****  St  *  *  $  *  *  *  ff  *  *  *  *  *  #  *  *  ft  ft  ft  ft  &  ft  ft  ft  ft  ff  ft  ft  ft  ft  ft  ft  ft  ftftftftftftftftffftftftftft 


SUBROUTINE  TWO(RHO.RX,RY.PX2.PY2,DX.N  ST2) 

I*  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *****  *  *  *  fc  *  *  4 


******* 


IMPLICIT  REAL*8(A-H.O-Z) 

DIMENSION  RX(N),RY(N),PX2(N).PY2(N) 
STWO  1=0.0 
STW02  =  0.0 
DO  20  1=  l.N-1 

TWOI  =  RX!l)*PX2(I) 

TWO  2  =  RY(IV:,PY2m 

stwoi  =  sT\voi  +  twoi 

STW02  =  STW02  +  TWO  2 
20  CONTINUE 

T21  =  (STWOl  +  RX(\)*PX2(N)  2.)*DX*RHO*2.0 
T22  =  (STW02  +  RY(N)*PY2(N)  2.)*DX*RHO  *2.0 
ST2  =  121  +  T22 
WRITE!  6.1 1 3  )ST2 
113  FORMAT(2X.  SI2  =  '.F 1 2.5) 

DU  I  VT»  'St  ft  ft  ft  ft  ft  ft  ft  ft  ft  ft  ft  ttffkftftft**1t*ftft1fttftfiftfHtflt*ftttftf 


DX*RHO*2.0 


PRINT*.'* 

RETURN 

END 


**********  St  :jt!!tSiiMt!f::|t!tt  St  StStStstStSt***  St  *******  **$**' 


******  **.**. St  ******  *******  ********  ********************************************* 

SUBROUTINE  THREE! N.RX.RY.PX1,PY1.PX2,PY2,TM,ST3,ST4,ST7, 

1  STS.ST11.ST12.ST16.sT1S) 


*  *  *  *  *  *  *  *  -is  *  *  *  *  *  *  *  *  *  *  *  v  *  $  *  *  *  : 


*  *  *  *  sjs  *  *  *  *•*  *  *  *  *  *  *  ******  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  ***** 


IMPLICIT  REAL*8(A-H.O-Z) 


DIMENSION  PX1(N).P\  1(NLPX2(N),PY2(N).RX(N).RY(N) 
ST3  =  (RX(N)*PX1(N)  +  RY(N)  PYllN )>*T\T-2. 

ST4  =  RX(N)'*PX2IN)  +  RY(N  >*PY2(N»*TM*2. 

CT7  _  B  V/  V  WDV 1/ V  .  D  V  1*PV  I  (V 


ST12  =  !PX[(N)-PY2(N)-PYl(NfPX2(N)):!TM 
ST  16  =  ( RX(N‘  >*PX  1(  N )  +  RY(  N  r  PY1<  N))*TM 
ST1 S  =  ( R  X  f  N  r  PX  2(N )  +  R  Y(N  r  P  Y2(  N )  r  T  M 
WRIT B6.1 1-4)ST3.ST4.ST7.ST8.ST1  US T12.ST16 
114  FORMA  I72X.ST3  =  .F12.5.2X.  ST4=  '.F12.5.2X.ST7=  .F12.5.2X,'ST8  = 

1  F12.5.2X.  ST1 1  =  '.F12.5.2X.  ST12  =  '.FI  2.5.2X.'ST16  =  '.F12.5.2X) 

WRITE!  6.1 1S)ST1S 
1  IS  FORMAT! 2X.  ST1S=  '.F12.5) 

p  ^  "J  *  *‘  '  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  f?  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  ' 

RETURN 

END 


^  *****  * 

C 


* *  * * *************** *  * *  *  * * *  *  *  * *  ****** * *  * *  *  *  *  *  * *  *  *  *  * *  * *  *  *  * * * *  *  *  *  * * *  * * *  *  * 

SUBROUTINE  FOU  RfN.RHO.DX.RX.RY.PXl.PYUPX2.PY2.ST5.ST9) 

* *  *  *  *  *  *  *  *  *  *  *  *  *****  * *  * *  *  * * *  *  * * *  *  * * *  * * * *  * * *  * * * * * * * * * * *  * *  * * * * * * * * *  * *  * *  *  * * 


IMPLICIT  REAL-S(A-H.O-Z) 
DIMENSION  RX(N).RY(N).PX1(N; 
SFOU1  =0.0 
SFOL2  =  0.0 
S  FOU  3  =  0.0 
SFOL  4  =  0.0 
DO  30  1=  UN-1 

FOU  1  =  R.X(I  )*PYl(h 
FOU2=  RY  IV-PXim 
SFOLT  =  SFOLT  -6  FOl  l 
SLOU2  =  SFOL  2-i-  LOU 2 
IOU3  =  RX<  I  )*PY2(l  > 
IOU4=RY(lV:PX2(I) 


1(N).PY1(N).PX2(N).PY2(N) 


on  no  nn  on 


SFOU3  =  SFOU3  4-  FOL‘3 
SF0U4=  SFOL4+  FOL'4 
CONTINUE 

T41  =  (SFOU 1  +  RX(N)*PYHN)  2.)*DX*RHO 
T42  =  S FOU  2  +■  RY(N)*PX  1  ( N  2.)*DX  *  R H 0 
T43=  SFOU3  +  RX(N)*PY2(N)  2>DX*RIIO 
T44  =  (SFOU4  +  RY(N)*PX2(N).  2.)*DX*RHO 
ST5  =  T41-T42 
ST9  =  T43-T44 


WRITE(6.1 15}ST5.ST9 

FORMAT(2X.'ST5  =  .F12.5.2X,ST9  =  ,  FI  2.5) 

PRINT"  ' #  *  *  *  *  *  *  *  *  *  #  *  *  *  *  *  *  *  *  $  *  *  *  *  *  *  *  *  *  *  *  ******  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  > 


PRINT*. 

RETURN 

END 


******  ***************************************************  ******..£  .j,  *  *  *  *  .j, 


SUBROUTINE  FIVE{N.RHO.DX.RXJRY.PX1.PYLPX2.PY2.ST6.STIO) 

******  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  .*  *  *  *  *  *  *  *  *  *  *  *  *  *  .*  *  *  *  * *  *  *  £  .*  #  *  *  *  * 


IMPLICIT  REAL*8(A-H,0-Z) 

DIMENSION  RX(N),RY(N),PX2(N),PY2(N),PX1(N),PY1(N) 

SFI V  1  =  (J.O 


SFIV2  =  0.0 
SFIV3  =  0.0 
SFIV4-0.0 
DO  40  1=  l.N-1 

FIV11-PX2(I)*PY1(I) 
FIV12=  PY2(IrPXli  I) 
FIV21  =  PXl(IrPY2i  I) 
FIV22«PYl]l)*PX2iT) 
SFIVl-SFIVI  +  FIVn 
SFIV2  =  SFIV2+  FIV12 
SFIV3  =  SFIV3  +  FIV21 
SFIV4=SFIV4+FIV22 
CONTINUE 


TM  =  (SF  VI  +  PX2(N)*PY1(N)  2.)*DX*R  IO 
Ti2=  SF  V 2  +  PY2(N f'PX  1(N j  2.}*DX*RHO 
T53  =  (SF  V 3  +  PXl(N)*PY2(N)  2.)*DX*RHO 
T54  =  (SF  V4+  PY1(N)*PX2(N)/2.)*DX*R  40 
SI6  =  T5I-T^2 
ST  10=  T53-T54 
WRITE(6.I16)ST6.ST10 

FORMAT(2X,  SI6  =  '.F12.5.2X  ,ST10  =  '  F12.5) 

PRINT*  **  *  *  *  *  *  *  *  *  *  *  *  .*  *  *  #  *  #  *  *  *  *  .*  *  *  *  *fi  «»:  $  *  *  *•  *  *  *  *  *  *  *  *  *  *  A 


PRINT-. 

RETURN 

END 


S  *  *  *  *  *  *  *  *  * ' 


*  *  *  *  *  **************************  *  *  *  ************ *  **********  * ***************** 


SUBROUTINE  SIX(N\RHO.DX.TM,RX.RY.PXl,PYl.PX2.PY2.ST19. 
ST20.ST21.ST22) 


*  vi  i  vj  i  _  r 

*i  *  *  *  .■{:  *  *  *  *  *  *  *  *  *  *  *  ******  *  *  *  *  *  *****  *  *  *  *  *  *  *  *  *  *  *  *  *  ********  *  *  *  *  *  *  *  *  *  *****  *  *  *  *  *  *  * 


IMPLICIT  REA L*8{ A-H .O-Z) 
DIMENSION  RX(N).RY(N).PX1(N). 


DIMENSION  RX(X),RY(N).PX1 
SS!XU  =  0.0 
SSIX12  =  0.0 
SSIX21  =0.0 
SSIX22  =  0.0 
SSIX31  =0.0 
SSIX32  =  0.0 
SSIX41  =0.0 
SSIX42  =  0.0 
DO  50  1=  l.N-1 

SIX1 1  =  PX2(I)*PYI(! 
SIX 1 2  =  PY2(  IrPXlU 
SIX21  =  RXdUPXKI) 


PY1(N).PX2(N).PY2(N) 


SIX11  =  PX2(I)*PY1(I) 
SIX  1 2  =  PY2(  IrPXlh) 
SIX21  =  RX(  I  )-PX  1(1 ) 
SIX22=  RY{I  rPYl(I) 
SIX31  =  PXi(i)::PY2(I) 
SIX32=PY1  I)*PX2(l) 


SIX41  =  RX(I)*PX2(  L) 
SIX42=RY<T)*PY2<T> 
SSIXU  =  SSIX11  +  SIX11 
SSIX12=  SSIX12  +  SIX12 
SSIX21  =  SSIX21  +SIX21 
SSIX22=SSIX22  +  SIX22 
SSIX31  =  SSIX31  +  SIX31 
SSIX32=  SSIX32  +  SIX32 
SSIX41  =  SSIX41  +  SIX41 
SSIX42  =  SSIX42  +■  SSX42 
50  CONTINUE 

T61  =  (SSIX1 1  +  PX2(N)*PY1(N) 


T61  ==  ( SSI X 1 1  +  PX2(N)*PY1(N)  2.)*DX*RHO 
T(>2  =  ( SS I X 1 2  •+-  PY2(N')*PXHN)  2>DX*RHO 
163  =  (SSIX21  +  RX(N)*PX1(N)  2.)t:DX*RHO 
T64=  SS1X22+  RY<N  *PY1(N  2>DX*RHO 
T65  -  ( SS  I X  3 1  +  PX I  ( N  f  P  Y2(  N)  2.  )*  DX  *  R  MO 


ST  19=  T61-T62 
ST20-T63  +  T64 
ST21  =  T65-T66 
ST22  =  T67  +  T6S 

WRITE!  6. 1 1 7jST19.ST20.ST21.ST22 


WRITE!  6. 1 1 7jST  19.ST20.ST2 1.ST22 
FORMAT(2X.  ST19=  '.f-T2.5,2X.'St20=  .F12.5.2X,  ST21  =  '.F12.5. 

12X.  ST22=  .F12.5) 

pD|^-  T  *  '  ❖  #  *  *  #  *  * &  *  «  *  *  *  #  *  *  *  *  *  *  *  *  *  *  #  #  *  *  *  *  *  *  *  *  sfs  *  *  ft  ft  *  *  *  *  *  *  *  *  *  *  ft  ft  *  ft  *  ft  ft  ' 

RETURN 

END 


PROGRAM  OF  DOUBLE  LINK  FLEXIBLE  BOOM  IN  3  DIMENSIONAL 
MOTION  WITH  TWO  MODES 


*  *  *  *  *  *  * ft  *  ft  ft  ft  ft  ft  ft  ft  ft  ft  ft  ft  ft  ft  ft  ft  ft  ft  ft  .-It  ft  ft  ft  ft  ft  ft  ft  ft  ft  ft  ft  ft  ft  ft  ft  ft  ft  ft  ft  ft  ft  ft  ft  ft  ft  ft  ft  ft  ft  ft  ft  ft  ft  ft  ft  ft  ft  ft  ft  ft  ft  ft  ft  ft  ft  ft  ft 

*  SIMULATIONS  OF  DOUBLE  LINK  FLEXIBLE  BOOM 

*  IN  3  DIMENSIONAL  MOTION  WITH  2  MODES 

j. 

ft  ft  ft  ft  ft  ft  ft  ft  ft  ft  ft  ft  ft  ft  ft  ft  ft  ft  ft  ft  ft  ft  ft  ft  ft  ft  ft  ft  ft  ft  ft  ft  ft  ft  ft  ft  ft  ft  ft  ft  ft  ft  ft  ft  ft  ft  ft  ft  ft  ft  ft  ft  ft  ft  ft  ft  ft  ft  ft  ft.  ft  ft  ft  ft  ft  ft  ft  ft  ft  ft  ft  ft  ft  ft  ft 

TITLE  SOLUTION  OF  SIMULTANEOUS  DIFFERENTIAL  EQUATIONS  FOR 
TITLE  DOUBLE  LINK  FLEXIBLE  BOOM  IN  3  DIMENSIONAL  MOTION 


TITLE  DOUBLE  LINK  FLEXIBLE  BOOM  IN  3  DIMENSION. 
T  IT  LE  WITH  TWO  MODES 


*  THIS  PROGRAM  USES  2  MODE  SHAPES  FOR  3D  AL.  AND  COM. 
FIXED  IER.  I  PVT.  N,  I 

CONST  DX  =  24..DELT  =  .04  .RHO  =  3.069UE-04.N  =  13 
PA  RAM  Cl  =  .4  ,C2  =  10.  ,C3  =  24.65 
INCON  X0=  3.1416 


INITIAL 


TM  =  9.7P6E-02 
RIBZZ=  145.39 
RIRZZ=  7.764 


'HIE  DAMPING  FORCE  WILL  BE  CHANGED 
BY  VARYING  THE  DAMPING  COL  I  I  I  ENT  ZETA 


OQ 


ZETA  =  0.0 
ZETA  =  0.002 
ZETA  =  0.005 


*  THIS  VALLES  ARE  FOR  ALUMINUM  ALLOY 


RX1L=  1.381 823E-02 
RZI L  =  -6. 144I72E-03 
RY2L=  1.527773E-02 


*  THIS  VALUES  ARE  FOR  COMPOSITE  MATERIAL 


RX1L  =  -1.395S93E-02 
RZ1L  =  -5.647S70E-03 
RY2L=  1.527779E-02 


DIMENSION  SIZE  SHOULD  BE  EXPRESSED  BY  NUMBER 
INSTEAD  OF  CHARACTER 


AKKA  i  lrVl(3),D(3) 


REMOVE  ASTRICKS  IN  DATA  STATEMENT  FOR 
EXECUTION  OF  EACH  MATERIAL 


DATA  RX  8.2085.16.417.24.6255.32-834.41.0425.49.251,57.4595, 

#  44.0388, 30.6 1S2, 17. 1976, 3. 777, -9. 6436, -23.(k>42 

DATA  RY  22.5526.45.1052.67.6578.90.2104,112.764,135.316,157.869, 

#  177.7653,197.662,217.559,237.456,257.353,277.25/ 


=  =  =  =  =  DATA  FOR  ALUMINUM  ALLOY  =  =  =  = 


DATA  OMG1  3.384515/, OMG2/3.568821/ 


DATA  PY 


-2.3202S8E-02. 

8.960620E-02. 

1.944046E-01, 

3.32801 7E-01. 

5.0O0194E-01. 

6.913154E-01, 

9.019S64E-01. 

1 . 1 990 10E  ■+•  00, 

1. 51 6552  E  +  00. 

1.S50174E +  00, 

2.195575E  +  00, -2.548629.-2.905414/ 


DATA  PX2  2.632094E-02  . 

=  1 .024 16  IE-0 1 . 

=  2.239926E-01, 

=  3.S676S6E-01. 

=  5.S64S63E-01. 

n  8.189327E-01. 

=  1.07994SE  +  00, 

=  1.331760E  +  00. 

p  1.601 282E  + 00. 

=  I.884697E  +  0O, 

=  2. 1 7S297E  +  00.2.4785 1 8.2.78 1 967; 


DATA  PZ2 


-9.574829E-03. 
-3.72661  KjE-02 
-S.1511 16E-02 
-1.401517E-01 
-2.1343SOE-01 
-2.98035  IE-01 


•:tucccccooccc;aoaococooc:o  ccocqoccocco 


* 

* 

* 

* 

* 

* 

* 

* 


*  -3.930303  E-01. 

#  -2.2317S7E-OI. 

3  -4.13S60SE-02. 

3  1.497749E-01. 

si  3.478057E-01.5.503027fe-01.7.549778E-01 


=  =  DATA  FOR  COMPOSITE  MATERIAL  =  =  =  =  = 


DATA  OMG1  4.287485  ,OMG2  4.461334 
DATA  PY1  2.38435<>E-02. 

*  9.20775  IE-02. 

9  1.997595E-01. 


3.419566E-01, 

5.137555E-01, 


9  7. 10279SE-01, 

=  9.266955E-01. 

=  1.21S125E  +  00, 

3  1. 530611 E +  00. 

=  1.S59596E  +  00. 

3  2.200665E  +  00.2.549586,2.902334/ 

DATA  PX2  2.632090E-02  , 

#  1.024162E-01, 

9  2.239929E-01, 

9  3.867693E-01, 

«  5.S64875E-01, 

3  8.1S9344E-01. 

3  1.07995  IE  +  00, 

3  1.33 1 764E  +  00, 

3  1.6012S6E  +  00, 

3  1.8S4702E  +  0O. 

3  2. 178303E  +  00, 2.478525 J.  78 1976/ 

DATA  PZ2  -9.575 142t-03. 

3  -3.726674E-02, 

3  -8.151225E-02, 

3  -1.407533E-01. 

3  -2.134401E-01, 

3  -2.980377E-01. 

3  -3.930335E-01, 

3  -2.23181  IE-01. 

3  -4. 138775E-02, 

=  1.497741  E-01. 

3  3.47S057E-01.5.503036E-01.7.549797E-01 

CALL  CONST  (N.DX.RHO.TM.RX.PY1  .PX2.PZ2.SST1.SST2.SST3.SST4.  ... 
SST5.SST6.SSt7.SST8) 

DERIVATIVE 

NOSORT 

*  TAO=_10_.-10.*STEP(6.S5) 

TACM  =  Cl  "‘TIME* *2 

TA02-C2*X1-  STEP(C3  )) 

TAO  =  SWITCH(TI  ME.  LE.T  I  ,TAO  1  .TA02) 

*  TAO  =  0. 

A(  1.1)=  RIBZZ+  RIRZZ  +  X2**2  +  X3**2*(  l.-2.*SST3-2.*SST4)  ... 

+  2.*X3*(SST1  +  SST21 
A(  1 .2)  =  SST7  +  SST8  +  ( SSI  5  +  SST6)*X3 
A(i.3)=-X2*(SST5+SST6) 

A(2.1)  =  AU.2) 

A(  2.2)  =1.0  ' 

A(2.3)  =  0.0 
A  3.1  =  All. 3) 

A(3.2)  =  O,0 
A(3.3)  =  1.0 

B(  1)=  TAO-2.*XlD*(X2D*X2  +  X3D*(X3  +  SST1  +  SST2-2.*X3*(SST3  +  SST4) ) ) 
B(2)  =  X 1  D**2*X2-2.*X1  D*X3D*(SSY5  +  SST6)-OMGl**2*X2  ... 
-2.*ZETA*OMG  1  *X2D 
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B(3)=  X1D**2*(X3-SST1-SST2  +  2.*X3*(SST3  +  SST4)  +  2.*X1D*X2D*(SST5 
+  SST6)  -OMG2**2*X3-2.*ZETA*O.MG2*X3D 
CALL  DGEFA  (A.  3.  3.  IPVT.  IER) 

IF  (IER.NE.O)  GO  TO  1 12 
CALL  DGESL  (A,  3.  3.  IPVT,  B.  0) 


X1D=INTGRL(X0.B( 
X2D=  INTGRL(0..B(2 
X3D=IXTGRL(O..Bi3 
XI  =  INTGRL!0.,X1D) 
X2=INTGRL(O..X2D 
. —  -(o.,x3r 


X3=INTGRL(0.,X3D) 

TH-X1 
THD=  X1D 
Q1  =  X2 
Q1D  =  X2D 
Q2  =  X3 
Q2D-X3D 
WX  =  PX2(N)*Q2 
WY=  PYIfNrOI 
WZ  =  PZ2(NJ*Q2 

WXYZ  =  WX^  +  WY**2+  WZ**2 
W-SQRKWXYZ) 

XSLOP=  RX1L*QI*57.2957 
YSLOP=  RY2L*Q2*57.2957 
ZSLOP=RZlL*Gl  *57.2957 

‘  WRITE!  6J  20)RI  BZZXSTJ.ST2.ST3.ST4.ST7.ST8, ST  1 1  ,ST1 2,ST16,ST  18, 
SY5.Sr9.ST10.STl9.Sf20,ST21,ST2i 
■120  FORMAT)  2X,  1  S( F 1 2.5.2X)) 

RETURN 

112  WRITE(6.il4)  TIME. IER 
114  FORMAT!  2X.0  IER  *  ,17) 

■  CALL  ENDJOB 


*  WRITE(6.1 16)(Q(IF1  ),IF1  =  l.N) 
*116  FORMAT)  2X.3F8.3.5X) 

PRINT  TH.THD,  Ql,Q2.\Vx.WY.WZ 

/-  -v  "r-  »  r*ix  A  t— v  t  -4  r, 


s^x^7?w^^ai^^5^-xsix»p-YSLop-zsix>p-TAo 

SAVE  0.025. TH,THD.Ql.Q2,TAO,WX,WY,WZ,W,XSLOP,YSLOP,ZSLOP 

*  PA  RAM  Cl  =  .8  ,C2=  20.  ,C3  =  14.0 

*  PA  RAM  Cl  =  1.6  ,C2  =  40,  ,C3  =  8.7 

GRAPH  (G 11  ,NJ  =  7.LO=  -S.00.DE  =  TEK618.SC=  8.0  )  TIME!  NT  =  5.LE  = 


GRAPH 


GRAPH 


SEC  j.TH(UN  =  RAD.RU  =  1X1=1) ... 

,TH(Xk-OMIT.PO-d  .RU  =  2. LI  =  2)  ... 

,TH(AX  =  OM IT.PO  =  0.  ,RU  =  3Xl  =  3) 

(G1.NT  =  ULO  =  0.,DE  =  YEK618,SC  =  .3  )  TIME(XI  =  5,LE  =  10..UN- 
SEC  )JHD(PO  =  .0.  UN  =  RAD  SEC  ,RU=  1,LI=  1) ... 
.THD<K>=.0.AX  =  OMIT  ,RU  =  2X1  =  2)  ... 


.  I  HO)  PO=  .0.AX-  OMIT  .RU  -  2. LI  = 
.THD  PO  =  .0,AX  =  OM  IT  .RU  =  3X1  = 
GRAPH  (G2. NT  =  8X6  =  -.80  .SC  =  .200,  t)E  = 
TIME! NT  -  5,LE=  10.. UN  =  SECh.Ql! 


EK618)  ... 

* _  * i  \* ^  r*  t 


TIME(NI  -  5.LE  =  10..UN  =  SEdT.QKUN  =  IN  ,RU  =  l.LI  *  1).  ... 
Ql(PO  =  0.,  AX  =  OVf  IT.RU  =  2X1  =  2).  ... 

Ql! PO  =  Oj,  AX  =  OMIT.RU  =  3.LI  =  3) 

/  C'  *!  X' f  —  a  r  r\ _ f\  aa  cc _  st\t\  nr _ Trt/zio\ 


GRAPH TG3NT  =  1L0  =  0.00  tSC  =  .500  .DE  =  TEK618)  ... 

TIMBNT  =  5.LE  =  10..UN  =  SEC  ).Q2(UN  =  IN  ,RU  =  l.LI  =  1) 
,Q2(PO  =  0„  AX  =  OMIT.RU  =  2. LI  =  2)  ... 

,02)  PO  =  0.  AX  =  OM IT.RU  =  3X1  =  3) 

GRAPH  (G5.DE  =  tEK618.LO  =  -0.0.SC  =  .60  .NT  =  8  )  TIME(NT  =  5,U: 
XE  =  10.).WX(UN=  IN'  ,RU  =  l.LI  =  1) ... 


m 


)  TIME(NI  =  5, UN  =  SEC'  ... 


)  TIME(NI  -  5,UN=  SEC’  ... 


AVX! PO  =  0..AX  =  OMIT.RU  =  2.LI  =  2) ... 

AVX)  PO  =  0.AX  =  OM  IT.RU  =  3, LI  =  3) 

GRAPH  (G6.DE  =  T£K618.LO=  -0.4. SC  =  .10  .  Nil  =  8  )  TIME(NI  =  5, UN  =  SEC'  ... 
XE  =  10.).WY(UN  =  TV  .RU=  IXl=  1) ... 

.WY(PO  =  0..AX  =  OMIT.RU  =  2. LI  =  2)  ... 

,WY)  PO  =  0.  AX  =  OM  IT.RU  =  3X1  =  3) 

GRAPH  (G7,DE  =  TEK618,LO  =  -.0.SC  =  .15  )  TIME(NT  =  5.UN=  SEC  ,LE=  10.) ... 


,WZ(l'N*  iN'.NI-S.RU-  1X1*1) ... 

\VZ  PO  =  O..AX  =  0\f  IT.RL  =  2. LI  =  2)  ... 

.WZ(PO  =  O.jAX  =  OM  IT.RL"  =  3.LI  =  3} 

GRAPH  (GS.DE  =  TfeK618.LO  =  O.O.SC  =  .55.NT  =  8.RU  =  I,2)TIME(NI  =  5.UN  = 
,LE  =  10.).W(U\=  TN'.RU  =  1X1=  I)... 

.W(PO  =  O..AX  =  OMlT.RU  =  2. LI  =  2) ... 

.W  PO  =  0..AX  =  OMlT.RU  =  3X1  =  3) 

)  TIME(\I  =  5X’N  * 


GRAPH  (G9,DE  =  TEK61S.LO  =  -  ASC  =  .1,NT  =  8.RU  =  I 
.I.£  =  10.).XSLOP(UN  =  DfeG  .RL  =  1X1=1).  ... 
XSLOP(PO  =  0.  .  “■  “ 


_  .  .AX -OMlT.RU -2X1  *2) 

,XSLOP(PO  =  0..AX  =  OV1IT.RU  =  3,L  I  =  i) 
GRAPH  (G4.DE  =  TEK618.LO=-1.2.SC=.3.NI  =  S,RU  =  1 
XL  =  10.). YSLOPI UN -  DEG  .RU-  1X1=1),  ... 
YSLOP(PO  =  O..AX*  OMIT.Rl'  =  2.LI  =  2).  ... 
YSLOP  (PO-O..AX-  OMlT.RU  =  3.LI  =  3V 
GRAPH  (GZ.DE  =  TEK618.LO*  -.40.SC  =  1,NI  =  8.RU  =  1 
XE  =  lOAZSLOPCUN  =  DE-G.RU  —1X1=1),  ... 
ZSLOP(PO  =  O..AX  =  OMIT.RL  =  2X1  =  2),  ... 

ZS  LOP  PO  =  0. .AX  =  OM  IT.RL  =  3, LI  -  3) 

Gil)  APPLIED  TORQUE 
G12)  ANGULAR  DISPLACEMENT 
Gl)  ANGULAR  VELOCITY 
G2  GENERALIZED  DISPLACEMENT  Qt 
G3)  GENERALIZED  DISPLACEMENT  02 


SEC' 

SEC'  . 

)  TIME(NT  =  5.UN  =  SEC' 

)  TIME(NI  —  5, UN  =  'SEC' 


LABEL 

LABEL 

LABEL 

LABEL 

LABEL 


“LABEL  <G4)  GENERALIZED  DISPLACEMENT  Q3 
LABEL  (G5)  DISPLACEMENT  IN  X-PiRECRJON 


LABEL 
LABEL 
LABEL 
LABEL 
LABEL 
LABEL 
END 
STOP 

FORTRAN 


DISPLACEMENT  IN  Y-Dl RECTION 
DISPLACEMENT  IN  Z-DIRECTION 
MAGNITUDE  OF  DEFLECTION  AT  TIP  POSITION 
X-SLOPE  AT  THE  TIP  POSITION 
.  Y-SLOPE  AT  THE  TIP  POSITION 
)  Z-SLOPE  AT  THE  TIP  POSITION 


C 

C**** 


!:  *  *  *  *  *  ft  ft  *  ft  *  ft  ft  *  *  *  *  *  ft  *  ft  *  *  *  .'I*,  ft  ft  ft  ft  ft  ft  *  *  *  ft  ft  ft  .ft  ft  ft  ft  ft  ft  ft  ft  ft  ft  ft  ft  ft  ft  ft  ft  ft  ft  ft  ft  ft  ft  ft  ft  ft  ft  ft  ft  ft  ft  ft  ft  ft  ft  ft  ft 


SUBROUTINE  CONST  (N.DX.RHO.TM.RX.PY1.PX2.PZ2.SST1.SST2. 
SST3.SST4.SST5.SST6.SST7.SS  roS 


£  *  ft  *  *  ft  *  *  *  *  *  *  *  «  *  ❖  *  *  *  *  *  *  *  *  «  *  *  ft  *  ft  *  *  *  *  *  *  *  *  *  *  *  *  *  «  *  *  *  *  ft  ft  *  «  ft  *  *  *  «  *  *  *  *  *  *  *  *  *  ❖  *  *  «  *  *  «  *  *  «  ft 

c 


10 


20 


I M  P  L I C I T  RE  A  L  *8(A- H  .O-Z) 

DIM EN SION  RX(  1 3 ).PY1(  1 3 J.PX2(  1 3),PZ2(  1 3) 

ST  1  =0.0 
S  13  =  0.0 
STS  *0.0 
ST7-0.0 
DO  10  I  -  l.N-1 

Tl  =  Rtf(I>«PX2(I) 

T3-PZ2(I)**2 
T5  =  PX2II  )*PYI(I) 

T7-  RXlI)*PYl(I) 

STI-ST1  +  T1 
ST  3  =  ST3  +  T3 
ST5  =  ST5  +  T5 
ST7-ST7  +  T7 
CONTINUE 

PRINT*  ST1.ST3.ST5.ST7 

SST1  =  ($T1  +  RX(N)*PX2(N)  2.0)*RHO*DX 

SST3  =  (ST 3  +  PZ2(N)**2  2.0)*RI  IO*DX 

SST5  =  ST5  +  PX2(N  )*PY  1  ( N )  2.0 1* R HO* DX 

SST7  =  (ST7^  RX(N)*PY1(N)  2.U)4RHO*DX 

SST2-T\I*RX(N)*PX2(N) 

SST4-TM*PZ2(N)**2 

SST6-TM*PX2IN)*PY1(N) 

SST8-TM*RX(N)*PYKN) 

WRITE!  6.20JSST1.SST '2.SST3.SST4.SST5.SST6.SST7.SST8 
FORMAT!  2X.S(I12.4.2X)) 

RETURN 
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END 


3.  PROGRAM  OF  DOUBLE  LINK  FLEXIBLE  BOOM  IN  3  DIMENSIONAL 
MOTION  WITH  THREE  MODES 


ft*^**#****#*^**^^*#****************#************#********  $#*##***  $***$*$*$#* 


SIMULATIOMS  OF  DOUBLE  LINK  FLEXIBLE  BOOM 
IN  3  DIMENS  ONAL  MOTION  WITH  3  MODES 


TITLE  SOLUTION  OF  SIMULTANEOUS  DIFFERENTIAL  EQUATIONS  FOR 
-  - LINK  pLEXIBLE  goOM  IN  3  DIMENSIONAL  MOTION 


TITLE  DOUBLE _ 

TITLE  WITH  THREE  MODES 


THIS  PROGRAM  USES  3  MODE  SHAPES  FOR  3D  ALUMINUM  ALLOY 


FIXED  IER,  I  PVT,  N.  I 

CONST  DX  =  24..dELT=.04  ,RHO  =  3.0690E-04,N=  13 

i  r\  i  \  (  «— >  «  «  *v  t  a  /-<  ^  i  /  r 


PA  RAM  Cl  =  .4  C2  =  10.  ,C3  =  24.65 

*INCON  X0=  1.570S 

❖ 


INITIAL 

TM  =  9.7176E-02 
ZETA  =  0.0 

*  ZETA  *  0.002 

*  ZETA  *0.005 

RIBZZ*  145.39 
RIRZZ*  7.646 
RX1L*  1.381  S23E-02 
RZ1L  =  -6. 144172E-03 
RY2L*  1.527773E-02 
RY3L  *  -3. 149467E-02 


*  DIMENSION  SIZE  SHOULD  BE  EXPRESSED  BY  NUMBER  INSTEAD  OF  CHARACTER 

* 

D  DIM ENSION  A(4,4),RX(  1 3),PY1(  1 3),PX2(  1 3),PZ2(  1 3),PX3(  1 3).PZ3(  1 3 ) 

ARRAY  IPVT(4),B(4) 

D  DATA  OMG1,  3.384515  ,OMG2, 3.568821  .OMC 3, 1.796723E  + 01 


D 

D 

D 

D 

D 

D 

D 

D 

D 

D 

D 


DATA  PY1  -2.320288E-02, 

«  -8.960620E-02, 

n  -1.944046E-01, 

=  -3.3280 1 7E-0 1 , 

4  -5.000 194E-01, 

=  -6.913154E-01, 

a  -9.01 9864E-0 1 . 

s  -1.199010E  +  00, 

a  -1.516552E  +  00, 

a  -1.850174E  +  O0. 

a  -2. 195575E  +  00,-2.548629,-2.905414; 


D 

D 

D 

D 

D 

D 


DATA  PX2 


2.632094E-02  , 
1.024161  E-01. 
2.239926E-01, 
3.867686E-01, 

5. 864863  E-01, 

8. 1S9327E-01, 
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■co  ccccccocccc  ccccccocccc  ccccccccooo  cacao 


# 

d 


1.079948E  +  00. 

I.331760E  +  00, 

1.6012S2E  +  00. 

1.S84697E  +  00. 

2. 1 78297E  +  O0.2.4785 1 8,2.78 1 967- 


DATA  PZ2; 


-9.574829E-03, 

-3.726606E-02, 

-8.151 1 16E-02. 

-1.401517E-01, 

-2.1 343S0E-0 1 , 

-2.980351E-01. 

-3.930303E-01, 

-2.231787E-01, 

-4. 13S60SE-02, 

1.497749E-01. 

3.478057E-01 ,5.503027fc-01 .7.549778E-01 , 


DATA  PX3'  1. 33889 IE-01, 

#  4.792555E-01, 

d  9.51 8909  E-01, 

d  1.46SS96E-01, 

*  l. 94988 1E-0O, 

d  2.318285E-00. 

d  2.5025S3E-00, 

tt  2.45904SE-00. 

d  2.226467E  +  ()(), 

d  1.83962  IE +  00. 

d  1. 335901 E  + 00, 7.5430 12, 1. 34 1246/ 


DATA  PZ3 


-4.S90S16E-02, 

-1.747S80E-01, 

-3.46990 1  E-01, 

-5.3534 17E-01, 

-7. 105826E-01, 

-8.44S452E-01, 

-9. 120993  E-01, 

-9.414692E-01. 

-1.098786E  +  00, 

-1.3599ISE  +  00, 

- 1 .699S70E  +  00,-2.092339,-2.5 108 1 7/ 


DATA  RX  8.2085,16.417.24.6255.32.834,41.0425.49.251,57.4595, 


=  44.038S.3O.6182. 17. 1976.3.777.-9.6436.-23.0642 

DATA  RY  22.5526.45.1052.67.6578.90.2104.112.764.135.316.157,869. 
*  177.7653,197.662.217.559.237.456.257.353,277.25, 


CALL  CONST  (N.DX.RH0.TM.RX.PY1,PX2,PZ2,PX3,PZ3.  ... 
SST1.SST11.SST2.  ... 

SST21,SST3,SST311SST4,SST4 1  .SST5.SST5 1  .SST6.SST6 1 
,SST7,S$T8.S$T9,SSYlO.S$T12,SSTl3) 


tCH^TI  YIE.lE.Tl.TA01.TA02) 


DERIVATIVE 
NOSORT 

*  TAO=J0_.-10.*STEP(6.85) 

TAOl  =  Cl  *TIY1E**2 
TA02-C2*11.-  STEP1C3  JL 
TAO=  SWITCI 
16  TAO  =  0. 

A(  1,1 ) -  RIBZZ  +  RI RZZ  +  X2**2  +  X3**2  +  X4**2  +  2.*X3*(SST1  +  SST2) 
+  2.*X4*(SST1 1  +  SST21>2.*X3**2*SST3-2.*X4**2*SST31  ... 
-4.*X3*X4*SST9-2.;TX3~*2*SST4  +  X4**2*SST41 )  ... 
-4*X3*X4*SST10 

A(  l  .2)  =  SST7  +  SST8  +  (SST5  +  SST6)*X3  + (SST5 1  +  SST6 1  )*X4 
A(  1 .3)=  -X2*(  SST5  +  SST6) 

A  1.4  =-X2*(SST51  +  SST61) 


A(2.3)  =  0.0 
A  2,4)  =  0.0 
A(3,l  )= A(  1,3) 

A  3  2)  =  0.0 
At  3,3)  =  t.O 
A  3.4)  =  0.0 
A(4,1)=A(1,4) 

A{4.2)  =  0. 

A  4.3  =0. 

A(4.4)=  1. 

B(l)  =  TAO*2.*XlD*X2D*X2-2.*XlD*X3D*(X3  +  SSTl  +  SST2-2.*(... 

X3*(SST3  +  SST4}  +  X4*{SST9  +  SST10)))-2.*X  1  D*X4D*{X4  +  ... 

~  - - -•*J-i"'9+SST10)  +  X4<!(SST31 +  SST41))) 

SST5  +  SST6)-2.*X  1  D*X4D*(SST5 1  +  SST6 1 )  ... 
_ >MG1*X2D 

B(3)=  X1D**2*(X3-SST1-SST2  +  2.*X3*(SST3  +  SST4)+  2.*X4*(SST9  +  SST10)  ... 

\  +  2.*X1D*X2D*(SST5  +  SST6»-OMG2**2*X3-2.*ZETA*OMG2*X3D 
B(4)=  X1D**2*(X4-SST1  1-SST21  +■  2.*X3*(SST9  +  SST10)  +  2.*X4*(SST3 1  +  ... 
SST41 ))+ 2.*X1D*X2D*(SST51  +  SST61)-OMG3**2*X4  ... 
-2.*ZETA*OMG3*X4D 
CALL  DGEFA  (A.  4.  4.  IPVT,  1ER) 

IF  (IER.NE.O)  GO  TO  112 
CALL  DGESL  (A.  4.  4.  IPVT,  B,  0) 

X1D  =  INTGRL(0.,B(  1" 

X2D=  INTGRL(0..B(2 
X3D=INTGRL  0..B  3 
X4D-INTGRL  (0..BM 
XI  =  1NTGRL(0..X1D) 

X2=INTGRL(0..X2D 
X3=I\TGRL  0.,X3D 
X4=INTGRL  0.,X4D 
TH-X1 
THD=  X1D 
~>1  =  X2 
>1D  =  X2D 
)2  =  X3 
)2D  =  X3D 
;3  =  X4 
)3D  =  X4D 
VX=PX2(N 


WV-PVK 


2  +  PX3(N)*Q3 


\VZ=  PZ2(N)*Q2  +  PZ3(N)*Q3 
■ . ~>  \VX  '*2  +  WY**2  +  w: 


Z**2 


1L*QI  *57.2957 
'2L*Q2  +  RY3L*Q3)*57.2957 
1 L  *01*57.2''''' 


WXYZ  ..  _ 

W  =  SO  RT(  WXYZ) 

XSLOP=RXlL'“ 

YSLOP= (RY 
ZSLOP=  RZ1  _ 

*  \VRITE(6.120)RIBZt..^  .  .  ... 

*  ST5.ST9.ST  1 0.ST  1 9.St20,ST2 1  ,ST2i 
*120  FORM AT(5X,  1 8( F 1 2.5, 2X )) 

RETURN 

112  WRITER!  14)  TIME.IER 
114  FORMAt(2Xj0  IER  =  ',I7) 

*  CALL  ENDJOg 

*  WRITE(6.1 16 
"1 16  FORMAT!  2X.j 


.  _  .2957 

ST1,ST2.ST3.ST4.ST7,ST8,STI 1.ST12.ST16.ST18. 


.  .,.IF1=  l.N) 

8,32X1 


PRINT  TH.THD.  Q  1.02. \VX.WY.WZ.W.XS LOP. YSLOP.ZSLOP.TAO 
CONTRL  £INTIM  =  40..DELPRT  =  .SO 

SAVE  ().025.TH.THD.Q1,Q2.TAO.WX,\VY.\VZ.W.XSLOP.YSLOP.ZSLOP 
END 

PA  RAM  Cl  =  .8  ,C2  =  20.  ,C3  =  14.0 
END 

PA  RAM  Cl  =  1.6  .C2  =  40.  .C3  =  8.7 

GRAPH  (G1  l.NI  =  7,LO  =  -8.00.DE  =  TEK618,SC  =  8.0  )  TIME(NT  =  5.LE  =  10. 


GRAPH 


LN  = 


SEC  ).TH(l. N  =  RAD,  RE  =  1,1.1=  1) 
,TH(AX  =  OM  IT.PO  =  <1  .RE  =  2,  LI  =  2) 


EN  = 
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,TH( AX  =  OM  I T.PO  =  0.  . RU  =  3.L  l  =  3 ) 

GRAPH  (GLXT  =  LLO  =  0.,DE  =  tEK6llSC  =  .3  >  TIME(M  =  5,LE=  10..UX- 
*SE(^)JHlJ(PO».6.X='RADSfeC.RU-l.Ll-T)... 
.THD(kO=.O.AX  =  OMIT.RU  =  2,LI  =  2) ... 

.THDCPO  =  .OAX  =  OMIT.RU  =  3X1  =  3j 
GRAPH  (G2. XI  =  S,L6  =  -.16  .SC -.040,  bE=  I  EK6I8)  ... 

TIME(M  =  5,LE=10..LN  =  'SEC^).Q1(LN=  r 
,Ql(PO  =  0.,  AX -OMIT.RU  -  2X1  =  2) ... 

.Q1(PO  =  0..  AX  —  OMIT.RU  —  3. LI  =  3) 
tGVNI  =  S.LO-0.00  .SC  =  .200  ,DE  = 


:-TVfRU-l.LI-l) 


GRAPH 


GRAPH 


(G3.XI  =  8.LO  — 0.00  .SC  =  .200  .DE  = 
,T I  \t  E( XI  =  5. LE  =  1 0.X  X  -  SEC').Q2( LX  = 
.Q2(PO  =  0.,  AX  =  O.VflT.RL  =  2.Lt  =  2)  ... 
.52(PO  =  0.  AX  =  OMIT.RL  =  3, LI  =  3) 
kj5.DE  =  tEK618.LO  =  -0.0.SC  =  .60  .XI  =  8 
.LE-  10.).WX(UX  =  IX'  .RU  =  1X1  =  1) .. 
,WX(PO  =  0..AX  =  OMIT.RL  =  2X1  =  2)  ... 
"X(PO  =  0.  **' 


TEK618) ... 

IX", RU  =  l.LI-1) ... 

)  TIME(XI  -  5.LX -  SEC' 


,\V 


.AX  =  OMIT.RL  =  3, LI  =  3) 

r*  «  *  y  I  A  •  ’  <  A  i  • 


GRAPH  (G6.DE—  f£K618.LO=  -0.4.SC  =  .10  .  XI  -  8  )  TIME(XI  =  5, LX  =  SEC' 
.LE-10.).WY(LX-'IX*  .RL=  l,Li=  1) ... 


,WY(PO  =  O..AX  =  OMIT.RU  =  2X1  =  2) 


TIME(XI  =  5, LX—  SEC'.LE-  10.) 


.VVYcPO  — 0..AX  =  OMIT.RL  =  3.LI  =  3 
GRAPH  (G7  DE  =  TEK618.LO=  -.0,SC  =  .15  ) 

,\VZ(LX-TX'.XI-8JRL-1,LI-1)  ... 

AVZ  PO  =  0..AX  =  OVilT.RL  -2.LI  =  2) 

.\VZ  PO- 0. AX  =  OMIT.RL  =  3, LI  =  3) 

GRAPH  (G8.DE  =  TfeK618.LO  =  6.0.SC  =  .55.XT  =  8,RL  =  1,2)TIME(XT  =  5, LX  =  SEC'  ... 
.LE  =  10.).W(LX=  IX'.RL  —  1X1  -  I)  ... 

,W(PO-  O. AX  =  OM ITJRL  =  5.LI  =  2) ... 

.WPO  =  0.  AX  =  OM  IT,kL  =  3.LI  =  3) 

GRAPH  (G9.DE  =  YEK618.L0  =  -  .4.S<?-.1,XI-8,RU  =  1  )  TIME(XI - 5X’X-  SEC'  ... 
,LE  =  10.).XSLOP( LX  =  DEG'.RL  =  1X1=  1),  ... 
kSLOP(PO  =  O..AX  =  OMrriRL  =  2.LI  =  2)  ... 

,  X  S  LOb  PO  =  0. .  AX  =  O  M  it.  RU  =  3X1  =  3) 

GRAPH  (G4.DE  =  TEK618.LO=-1.2.SC  =  .3.XT  =  8.RL=  1  )  TIME(AT  =  5.LX-  SEC'  ... 
.LE  =  10.).YSLOP(LX-  DE&.  RU  =  1X1=  1), ... 

VSLOP( PO  =  0. AX  =  OM  IT.Rt  =  2.Lf  =  2),  ... 

YSLOPLPO-O.AX-  OMIT.RL  =  3X1  =  3) 

GRAPH  (GZ.DE  =  TEK618XO- -.40.SC  =  l.XI  =  S.RU  =  1  )  TIME(X1  =  5.UX  =  SEC'  ... 
.LE  =  10. ).ZSLOP(LX=  DECj'.RL-  1X1=  1),  ... 

ZSLOP(PO  =  O  Aa  =  OMIT.RL  =  2X1  =  2),  ... 

ZSLOP  PO-O.AX- OMIT.RL  =  3X1  =  3 
LABEL  (G 11)  APPLIED  TORQUE 
LABEL  G12)  AXGLLAR  DISPLACEMEXT 
LABEL  Gl)  AXGLLAR  VELOCITY 
LABEL  G2)  GEXERALIZED  DISPLACEMEXT 
LABEL  (G3)  GEXERALIZED  DISPLACEMEXT 
•'LABEL  (G4)  GEXERALIZED  DISPLACEMEXT 


LABEL 

LABEL 

LABEL 

LABEL 

LABEL 

LABEL 


G5)  DISPLACEMEXT  IX  X-DIRECRIOX 


G6)  DISPLACEMEXT  IX  Y-DIRECTIOX 
G7  DISPLACEMEXT  IX  Z-DIRECTIOX 
GS  MAGXTTUDE  OF  DEFLECTIOX  AT  TIP  POSITIOX 
G9)  X-SLOPE  AT  THE  TIP  POSITIOX 
(G4)  Y-SLOPE  AT  THE  TIP  POSITIOX 
LABEL  (GZ)  Z-SLOPE  AT  THE  TIP  POSITIOX 
EXD 
STOP 
FORTRAX 

SUBROUTIXE  COXST  (X,DX.RH0.TM.RX.PY1,PX2.PZ2.PX3.PZ3. 
1  SST1.SST11.SST2. 

1  SST21  .SST3.SST3 1  .SST4.SST4 !  .SST5,SST5 1  .SST6.SST61 . 

1  SST7SST8.3ST9.SST  1 0.SST  12.SSTI3) 

IMPLICIT  Rfi.AL *8(A-H.O-Z) 

DIMEXSIOX  RX(13),PY1(  13).PX2(  13),PZ2(  13),PX3(  13).PZ3(  13) 

sri  1  =  0.0 

ST3  =  0.0 
ST3 1  =  <>.0 
ST  5  =  0.0 
ST5 1  =  0.0 
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ST7  =  0.0 
DO  10  1=  UN-1 


■  Rtf(I)*PX2(n 
-RXn>*PX3(t) 


T5=  PX2(IrPYl(I) 
T51  =  PX3(I)*PVi(i) 
T7=RX(1>PY1(1 


T7=RX(Ij*PYl(I 
T9=  PZ2(  IJ*PZ3(  I ) 
T13  =  RX(1) 
ST1-ST1+T1 
STll-STll  +  Tll 
ST3  =  ST3  +  T3 
ST31-ST31  +  T31 
ST5=  ST5  +  T5 
ST51  =  ST51  +  T51 
ST7-ST7  +  T7 
ST9  =  ST9  +  T9 
ST13  =  ST13  +  T13 
CONTINUE 


SST3  =  (ST3  +  PZ 2(N)** 2  2 .  Oj  *  R  H  O*  D  X 
SST31  =(ST31  +  PZ3(NV**2:2.0>*RHO*DX 
SST5  =  (ST5  +  PX2(NV PY1(N)  2.0)*RHO*DX 
SST51  =  (ST51  +  PX3(NT PY  1(N )  2.0 U RHO*  DX 
SST7-(ST7+RX(N>PY1(N)‘2.0)*RHO*DX 
SST9  =  ST9  +  PZ2(N)*PZ3(N  >  2.0)*RHO*DX 
SST13  =  (ST13  +  RX(  N')*  PX3(N )/  2.0  )*  RHO*  DX 
SST2  =  T\l*RX(N')*PX2(N) 

SST21  =  TM*RX(N)*PX3(N) 

SST4  =  TVl  *  PZ2(N );  *2 
SST41  =  TM*PZ3(N}**2 
SST6=  T\l*PX2(NrPYl(N) 

SS"61  =  T  M  *  PX  3(  N  i*  P  Y 1  ( N ) 


SST21  =  T\1*RX(' 
SS  T4  =  T\T*PZ20 
SST41  =  TM*PZ3( 
SST6=  TM*PX2(. 


C  WRITE(  6.20  )SST  1  .SST2.SST 

C  20  FOR\lA"l(2X.S(F12.4.2X)) 
RETURN  ’ 

END 


SST4,SST5,SST6,SST7,SST8 
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